NON-COMMUTATIVE REPRESENTATIONS OF FAMILIES OF k 2 
COMMUTATIVE POLYNOMIALS IN 2k 2 COMMUTING VARIABLES 



HARRY DYM, J. W. HELTON, AND CALEB MEIER 

Abstract. Given a collection V = {p\(x\, . . . ,x 2 k 2 ), ■ ■ ■ ,Pk 2 i x i, ■ ■ ■ , x 2k 2 )} of 
k 2 commutative polynomials in 2k 2 variables, the objective is to find a condensed 
representation for these polynomials in terms of a single non-commutative polyno- 
mial p(X, Y) in two k x k matrix variables X and Y. Algorithms that will gener- 
ically determine whether the given family V has a non-commutative representation 
and that will produce such a representation if they exist are developed. These algo- 
rithms will determine a non-commutative representation for families V that admit a 
non-commutative representation in an open, dense subset of the vector space of non- 
commutative polynomials in two variables. 



1. Introduction 

This paper addresses a new type of problem concerning a condensed description of 
a collection of polynomials. 

1.1. Problem statement. The problem is to represent a family V of k 2 polynomials 
Pi, ....,Pfe 2 in 2/c 2 commuting variables xi, x 2 k 2 as an nc (non-commutative) poly- 
nomial p = p(X, Y) in two k x k matrices X and Y whose entries are the variables Xj 
(without repetition). For example, it is readily checked that if 

Pl(xi, . . . , X$) = x\ + X 2 X 3 + XiX 5 + X 2 Xj 
p 2 (xi, ■ ■ ■ , X 8 ) = XiX 2 + X 2 X A + XiX 6 + X 2 Xg 

p 3 (x l ,..., X S ) = XiX 3 + x 3 x A + x 3 x 5 + x 4 x 7 
p 4 (xi, ...,x 8 ) = x 2 x 3 + x\ + X 3 Xq + X A Xg, 

then 

( Pl P2 ) =X 2 + XY with X = ( Xl X A and Y = h X °) . 

\P3 Pa) \x 3 x a J \x 7 x 8 J 

The main objectives of this paper are to: 
(1) Present a number of conditions that a given set of polynomials 

p 1 (x 1 , . . .,x 2k2 ), . . . ,p k 2(x 1 , . . .,x 2k2 ) 
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in 2k 2 commuting variables x±, . . . , x 2 k 2 must satisfy in order for it to admit an 
nc representation p(X, Y). 

(2) Present a number of procedures for recovering such representations, when they 
exist. 

To formally describe the problem we shall say that a family V of k 2 polynomials 

Pi = Pi{xi, ■ ■ ■ ,X 2 kz), ...,Pk = Ptfixi, ■ ■ ■ ,X 2 tf), 

in 2k 2 commutative variables xi, • • • , x 2 k 2 , admits a nc representation if there exists 
a pair ofkxk matrices X and Y and an nc polynomial p in two nc variables such that 



/ 



(1.1) 



X 



X a {2) 



X a (k) 



\ 



(1.2) 



and 



(1.3) 



Y 



\ %a(k(k-l)+l) a; cr(fc(fc-l)+2) • 

/ ^o-(fc 2 + l) ^cr(fc 2 +2) 



\ ^cr(fc 2 +fc(fc-l)+l) X a (k2+k(k-l)+2) 
( P\(l) P\(2) 



p(X,Y) 



X a (k 2 ) ) 

■ ■ a;<r(fc 2 +fc) \ 

• • ^cr(2fc 2 ) J 

■ ■ P\(k) \ 



\ P\(k(k-1)+1) PA(fe(fe-l)+2) • • PX(k?) J 

where a is a permutation of the set of integers {1, . . . , 2k 2 } and A is a permutation of 
the set of integers {1, . . . , k 2 }. 

1 .2. Our Algorithms. The main contribution of this paper is to introduce a collection 
of algorithms for solving the nc polynomial representation problem. Since they are 
long, full descriptions are postponed to the body of the paper. However, we shall try to 
present their flavor in this subsection. 

The algorithms are based on the analysis of the patterns of one letter words, two 
letter words and some three letter words in the given family of polynomials. Thus, for 
example, if the one letter word 7x§ occurs in one of the polynomials, and the family 
admits an nc representation p(X, Y), then x 5 must be a diagonal entry of either X or 
Y and k — 1 of the other polynomials will contain either exactly one or two one letter 
words of the same degree with the coefficient 7. If there are no one letter words, then 
the analysis is more delicate; see Sections 7 and 7.2. 

If the diagonal variables are determined successfully, then subsequent algorithms 
serve to partition the remaining 2k 2 — 2k variables between X and Y and then to 
position them within these matrices. En route the k 2 polynomials are arranged in an 
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appropriate order in a k x k array. The final step is to obtain the nc polynomial; this is 
done by matching coefficients as in Example 1 . 1 below. 

Most families of polynomials containing k 2 polynomials in 2k 2 variables will not 
have nc representations. Either there will be no way of partitioning and positioning 
the variables that is consistent with the family or there will be no choice of coefficients 
that work. 

1.2.1. Examples. The next examples serve to illustrate some of the structure one sees 
in families of polynomials V that admit an nc representation and how it corresponds to 
diagonal determination, positioning and partitioning. 

Example 1.1. The family of polynomials 

Pi = 3x 2 x 4 + 3x 4 x 8 + x 2 x 3 + x 4 x 7 + 6x1X3 + 6x3X7 + X1X4 + x 3 x 8 
p 2 = 3x 2 x 6 + 3x 6 x 8 + xix 6 + x 5 x 8 + 6x1X5 + 6x5X7 + x 2 x 5 + x 6 x 7 
P3 = 3x 2 + 2xiX2 + 3x4X6 + X4X5 + X3X6 + 6xj + 6x3X5 
Pi = 3x 4 x 6 + 3xg + 2x 7 x 8 + x 3 x 6 + 6x3X5 + 6x7 + X4X5 
admits an nc representation. 

Discussion: If the given family of polynomials admits an nc representation p(X, Y), 
then it must admit at least one representation of the form 

(1 .4) p(X, Y) = aX 2 + bXY + cYX + dY 2 

for some choice of a, b,c,d G R, since pi, . . . ,p 4 are homogeneous of degree two. 
Moreover, as 

P3 = 3x 2 + 6x1 + 2xiX2 + • • • and p 4 = 3x 8 + 6x7 + 2x7X3 + • • • 

and there are no other one letter words in the family V, it is not hard to see (as we 
shall clarify later in more detail in §2) that x 2 and x 8 are diagonal entries in one of 
the matrices, say X, and correspondingly x x and X7 are diagonal entries in the other 
matrix, Y . Moreover, since x 2 and x 6 are in p 3 , whereas X7 and x 8 are in p 7 , it follows 
that if an nc representation exists, then, either 

(1.5) X=(* I), J ? ) and *X,Y)={* , 
or 

(1.6) X=(* I), K=(7 I) and «X,Y)=§ , 

and, in (1.4) we must have 

a = 3, (b + c) = 2 and d = 6. 

We shall assume that (1.5) holds; the other possibility may be treated similarly. 

The next step is to try to partition the remaining variables between X and Y. To- 
wards this end it is useful to note that if 

X=( X2 Xa ) then x 2 =( X2 + XaXb ') 

\%b x 8 J \ ' "/ 
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and hence (since we are assuming that (1.5) is in force) that p 3 must contain a term of 
the form 3x a x&. Comparison with the given polynomial p 3 leads to the conclusion that 
£4 and x e belong to X. Let us arbitrarily position x e as the 12 entry of X and x 4 as 
the 21 entry of X. Then 

X2 Xq\ , v 2 _ ( x\ + X§Xt± ^2^6 + %%X% 



X = r z )i and X 

\X4 X 8 

The 11 entry of X 2 provides no new information, but comparison of the 12 entry with 
the given polynomials leads to the conclusion that if the given family admits an nc 
representation, then p 2 must sit in the 12 position in p(X, Y) and hence 



p(X,Y) 



P3 P2 
Pi Pi 



Similarly, 



Y 





Xc\ _ 


^Y 2 = 1 




X 7 J " 





which upon comparison with the entries in p 2 leads to the conclusion that x c = x$. 
Therefore, Xd = X3. Comparison of 



3X 2 + bXY + cYX + 6Y 2 with 



P3 P2 
Pi Pi 



implies further that equality will prevail if and only b = c = 1 (i.e., a = 3, b = c = 1 
mdd = 6 in (1.4)). 

The example illustrates the strategy of first determining which variables occur on 
either the diagonal of X or of Y. In general, if the given family V admits an nc 
representation p(X, Y) of degree d and if 

p(X, Y) = aX n + bY n + --- with \a\ + \b\ > 

(and no other nonzero multiples of X n and Y n ) for some positive integer n > 2, then: 

(1) if b = (resp., a = 0), there will be exactly k one letter words of degree n with 
coefficient a (resp., b); 

(2) if ab 7^ and a ^ b, there will be exactly k one letter words of degree n with 
coefficient a and exactly k one letter words of degree n with coefficient b; 

(3) if a = b, there will be exactly 2k one letter words of degree n with coefficient 

a 

Example 1.1 fits into setting (2). 

So far we have focused on how we can use one letter words occurring in polyno- 
mials in V. A substantial part of this paper is also devoted to developing procedures 
for finding the diagonal variables that are based on patterns in two and (some) three 
letter words. The latter come into play if there are no one letter words to partition the 
diagonal variables between X and Y. 
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1.3. Effectiveness of our algorithms. Our algorithms depend upon the existence of 
appropriate patterns of one, two and some three letter words in the k 2 polynomials 
in the given family. We shall show that these algorithms are effective generically, 
i.e., they are effective on an open dense set of the set of polynomials that admit nc 
representations. 

1.3.1. General Results. Let W be the space of nc polynomials in two variables of 
degree d. We say that a subspace hi of W is of degree d if the maximum degree of the 
basis elements of hi is d. 

Theorem 1.2. Letpi, . . . , p^ be a family V of polynomials in 2k 2 commuting variables 
x\, . . . , %2k 2 of degree d > 3 and let U be a subspace ofW of degree d. Then there 
exists an open dense subset Sofhl for which the algorithms developed in this paper 
determine an nc representation p 6 S for V if and only ifV has an nc representation 
p G S. If such a representation exists, then these algorithms may be used to construct 
it. 

Proof. The proof is postponed until Section 7.3. □ 

Example 1.3, below, exhibits a family V with an nc representation p for which the 
algorithms do not work. However, perturbing p gives V for which they do. 

Example 1.3. Suppose that we are given the list of polynomials 

Pl = X2X4 + X4X8 + X 2 X 3 + X4X7 + X1X3 + X3X7 + X1X4 + x 3 x 8 
p 2 = X 2 X e + X 6 X 8 + XxX 6 + X 5 X 8 + XxX 5 + X5X7 + x 2 x 5 + x e x 7 

p 3 = x\ + 2xix 2 + x 4 x 6 + £4X5 + x 3 x 6 + x\ + X3X5 
p 4 = x 4 x 6 + x\ + 2x 7 x 8 + x 3 x 6 + £3X5 + x 2 + £ 4 £ 5 

Discussion If the given family admits an nc representation p(X, Y), then it is readily 
seen from the one letter words in the family V that it must be of the form (1.4) with 
a = d = 1 and that x\, x 2 , £7 and x 8 are diagonal variables. 

However, it is impossible to decide on the basis of one letter words which of these 
variables belong to X and which belong to Y. The most that we can say so far is that 
x\, x 2 and p 3 must lie in the same diagonal position, and hence 27, x 8 and p^ are in the 
other diagonal position. Thus, we may assume that x 2 is in the 11 position of X, x\ is 
in the 1 1 position of Y . But it is still not clear how to allocate x 7 and x 8 . 

It is readily seen that p(X, Y) = X 2 + XY + YX + Y 2 is an nc representation for 
the family of polynomials in Example 1.3. Thus, p(X, Y) is contained in the subspace 
hi of nc polynomials defined by 

U = {aX 2 + bXY + cYX + dY 2 : a,b,c,de R}, 



and, although our algorithms are not effective on the entire subspace hi, they do work 
on the (open dense) subset S of hi consisting of polynomials in hi for which a ^ d and 

2a^b. ~ □ 
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1.3.2. More detailed statements. Our main theorems on algorithm effectiveness are 
more detailed than Theorem 1.2. These theorems and a number of our algorithms 
depend in part on the coefficients of the terms in the commutative collapse p of an nc 
polynomial p(X, Y), which is defined as the commutative polynomial 

p(x,y) = p(xl,yl). 

In particular, if (p(i,j) is the sum of the coefficients of the terms in the nc polynomial 
p(X, Y) of degree i in X and degree j in Y, then (p(i,j) is the coefficient of x l yi in 
p(x,y). 

We shall also need the following more refined quantities: 

(f(i,j; X) (resp. <p(i,j; Y)) denotes the sum of the coefficients of the terms in the nc 
polynomial p(X, Y) of degree i in X and degree j in Y that end in X (resp. end in Y); 

<p(X; (resp. ip(Y; denotes the sum of the coefficients of the terms in the nc 
polynomial p(X, Y) of degree i in X and degree j in Y that begin with X (resp. begin 
withF). 

Thus, for example, if 

p(X, Y) = aX 2 YXY + bXYXYX + cYXYX 2 + dYX 3 Y, 

then 

<p{X;3,2) = a+b, <p(Y; 3, 2) = c+d, ip (3,2; X) = b+c and y?(3, 2; Y) = a+d 
Clearly 

The next theorem provides some insight into our one letter algorithms. There is 
an analogous result for two letter words: Theorem 7.14, which will be presented in 
Section 7.3 

Theorem 1.4. Letp\, . . . , p^ be a family V of polynomials in 2k 2 commuting variables 
Xi, . . . , x 2 fc2 and let W denote the set ofnc polynomials p(X, Y) of degree d > 1 such 
that there exists an integer t > 2 for which 

¥>(t,0)^0, <^(0,t)^0, cp{t,0)^cp{0,t) 

and that additionally satisfies one of the following properties: 

(1) ^(t,0)^(t-l,l) 

(2) t<p(p,t)^<p(l,t-l) 

(3) <p(t - 1, 1; Y) ^ and <p(t, 0) ^ <p(t - 1, 1; Y) 

(4) ip(Y; t - 1, 1) ^ anJ ^(t, 0) ^ </j(V; t - 1, 1) 

(5) (p{l,t-l;X) jt0and<p(0,t) + (p(l,t-l;X) 

(6) (p(X; 1, * - 1) ^ anJ ^(0, t) + tp{X\ 1, t - 1). 

r/zen one fetter algorithms stated in Section 4.7 determine that V admits an nc 
representation p(X, Y) in the class W if and only if V has a representation in this 
class. If such a representation exists, then these algorithms can be used to construct it. 

Proof. See Theorems 4. 14, 4. 1 5 and Remark 4.16. □ 
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The long list of caveats looks confining, but they are all strict inequality constraints 
and so hold generically. 

1.4. Uniqueness. The issue of uniqueness of an nc representation is of interest in its 
own right. We shall see is that while the representation p(X, Y) is highly non-unique, 
the arrangement of commutative variables Xj in the matrices X and Y is determined 
up to permutations, transpositions and interchanges of X and Y. 

1 .4. 1 . Polynomial identities and non-uniqueness of p. A basic theorem in the theory of 
rings with polynomial identities implies that if £ 2 fc denotes the set of all permutations 
of the set {1, ... , 2k} for each positive integer k, then the polynomial 

(1.7) q(X u . . . ,X 2 k) = ^ sgn{a)X a (i) • • • -Xgpfc) = 

for every choice of the k x k matrices Xi, . . . , X 2 k in R. kxk . Thus, if X and Y are 
arbitrary real k x k matrices and if 

(1.8) P (x,Y)= f g (x 1 ,...,x 2k ) 

with 

(1.9) X 3 = [X\Y] forj = l,...,2k, 

then p(X, Y) = 0; see [3] and [1] for additional information. 

Any other replacement of Xj in (1.7) by a polynomial in X and Y will also yield 
a polynomial p(X,Y) = 0. However, the choice in (1.9) will have nonzero coeffi- 
cients. In particular this means that if a given family V of polynomials admits an nc 
representation, then it admits infinitely many nc representations. 

If k = 2, for example, the nc polynomials 

(1.10) YXY 2 X + Y 2 X 2 Y + YX 2 YX + XY 2 X 2 + XYXY 2 + X 2 YXY 
and 

(1.11) Y 2 XYX + YX 2 Y 2 + XY 2 XY + YXYX 2 + XYX 2 Y + X 2 Y 2 X 

generate the same family regardless of how the commutative variables x±, • • • , x& are 
partitioned between X and Y and positioned. 

This stems from the fact that the difference between the nc polynomial in (1.10) and 
the nc polynomial in (1.1 1) is equal to the commutator 

[Y - X, (XY - FX) 2 ] = (Y — X)(XY - YX) 2 - (XY - YX) 2 (Y - X) = 0, 

since for 2x2 matrices 

X = ( Xl X2 ) and Y - 1 ' 



the polynomial (XY — YX) 2 has the special form 

(1.12) p(X,Y) = [X,Y] 2 = (XY - YX) 2 = ^ 
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This is well known by the experts in matrix identities, and it is easily verified by direct 
calculation that 

p(x) = x\x\ - X2{x 3 [2xqx 7 + (x 5 - x s ) 2 ] - (xi - Xi)x 7 (x b - x 8 )} 

+ X^\X\X-iX^ — £3X4X5 + X3X6 — X ± Xj + 2X1X4X7 — X4X7 + X3[—Xi + XijXs}, 

Thus, the polynomial p(X, Y) in (1.12) is an example of a homogeneous nc polynomial 
that produces a family V with some of the polynomials in V equal to zero and some 
not. We shall see later that the fact that the degrees of the polynomials in V are either 
4 or is consistent with Lemma 4.7. 

If p(X, Y) = (X+Y ) n for some positive integer n, then it is impossible to determine 
which variables belong to X and which variables belong to Y. 

The theorems presented later in the paper that validate our algorithms, e.g., Theorem 
1.4, have hypotheses that exclude cases like (1.10). 

1.4.2. Uniqueness ofX, Y. We just saw that an nc representation for V is highly non- 
unique, however, the pair X, Y in such representations is generically very tightly de- 
termined. This is indicated by the following theorem. 

Theorem 1.5. IfV is a family of polynomials pi, - ■ ■ , Pk 2 in the commutative variables 
xi, . . . , X2k 2 that admits two nc representations p(X, Y) and p(X, Y) that satisfy the 
conditions of Theorem 1.4, then there exists a permutation matrix U such that one of 
the following must hold: 



(1.13) 



1) 


X 


= n T xn, 


Y 


= U T YU, 


2) 


X 


= n T x T n, 


Y 


= Tl T Y T Tl, 


3) 


X 


= n T rn, 


Y 


= n T xn, 


4) 


X 


= n T r T n, 


Y 


= n T x T n. 



Proof. The proof is postponed until Section 4.4. □ 

We shall say that the pairs X, Y and X, Y are permutation equivalent if they are 
related by any of the four choices in (1.13). 

1.5. Motivation. The problem we study in the paper is undertaken primarily for its 
own sake, however, the original motivation arose from the observation that the running 
time for algebraic calculations on a large family of commutative polynomials V can be 
much longer than the corresponding calculation on a small family of nc polynomials 
representing V . Such calculations can be done using nc computer algebra, for example 
NCAlgebra or NCGB [2], which runs under Mathematica. 

As an example, consider computing Grobner Bases, a powerful but time consuming 
algebraic construction. The reader does not need to know anything about them to get 
the thrust of this example. We have a list P of nc polynomials and run an nc Grobner 
Basis algorithm on 

(1.14) P = {a m + m a + m m, aw + w w + w a , m am, 

m T aw, m T a T m, w T aw, w T a T m, w T a T w}. 
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Using NCGB on a Macbook Pro it finished in .007 seconds. Now substitute two by 
two matrices 



with commuting entries for the variables and run the ordinary Mathematica Grobner 
Basis Command. The run took 161 seconds in the most favorable monomial order that 
we tried. The corresponding 3x3 matrix substitution yielded a Grobner Basis compu- 
tation which did not finish in 1 hour. Calculation with higher order matrix substitutions 
would be prohibitive. 

The NC Grobner Basis and the commutative one contain different information. 
Namely, the NCGB determines membership in the two sided ideal Xp generated by 
P while the commutative GB obtained from "generic" n x n matrix substitution de- 
termines membership in the ideal generated by Xp + R n where R n is the ideal of all 
nc polynomials which vanish on the n x n matrices. We would assert that the NCGB 
contains very valuable information (possibly more than in the X P + R n case) and is 
readily obtained. In fact what brought us to the nc polynomial representation question 
was the reverse side of this. To speed up nc GB runs we tried symbolic matrix substitu- 
tions in the hope that the commutative GBs would go quickly and as n got bigger guide 
us toward an NCGB. This approach seems hopeless because of prohibitively long run 
times. 

In special circumstances nc representations could have a significant advantage for 
numerical computation. In particular, the numerics for solving the second order poly- 
nomial (in matrices) equation, called a Riccati equation, is highly developed. Con- 
sequently, it would be very useful to be able to replace a collection of conventional 
polynomial equations by an nc representation. 

Finally we mention that there is a burgeoning area devoted to extending real (and 
some complex) algebraic geometry to free algebras. Here one analyzes non-commutative 
polynomials with properties determined by substituting in square matrices of arbitrary 
size. See the recent references [4, 5, 6, 8, 9] and their extensive bibliographies. 

1.6. Computational Cost. The problem considered here can be attacked by "brute 
force" rather than by the methods developed in this paper. There are (2/c 2 )! arrange- 
ments of the variables sci, * • • , x 2 k 2 m X and Y and (k 2 )\ arrangements of the poly- 
nomials pi, ■ ■ ■ ,p k 2 in V. For a given arrangement a of the variables in X and Y, 
one obtains a matrix of commutative polynomials by forming a general nc polynomial 
p(X, Y) = J2 a 8 c a/3 m ai3(X, Y) of degree d as in (2.3) with undetermined coefficients 
c a/ 3 that are chosen to match the array determined by the arrangement A of k 2 polyno- 
mials, if possible. For each pair of arrangements a and A, one attempts to solve for the 
coefficients c a p to obtain an nc representation. We will refer to this approach as the 
Brute Force Method. Because there are [k 2 )\{2k 2 )\ possible systems, the cost of this 
approach is very high. Also, to rule out the existence of an nc representation this way 
it is necessary to check all of these cases and to verify that they fail. 

Much to the contrary, the procedures we introduce are likely to determine non- 
existence in the first few step. Even when me must run through all the steps, we 
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find that the implementation of the algorithm that we call Algorithm 2 requires on the 
order of 

d 2 3i +! 

10 (k 7 + 3dk 5 + d 3 k 3 ) + V 

3 

i=2 



d 2 3i + 1 



operations, which is much less than the 

(2P)!(A; 2 )! 3 
operations required by brute force; see §6 for details. 

2. One and two letter monomials in the k 2 polynomials: Determining 

THE DIAGONAL VARIABLES 

In this section we enumerate the one and two letter monomials that appear in the 
k x k array of commutative polynomials corresponding to the nc polynomial p(X, Y). 

2.1. Preliminary calculations. This subsection is devoted to notation and a couple of 
definitions that will be useful in the main developments. 

Let e±, ... ,ek denote the standard basis for IR fc and let E st denote the k x k matrix 
with a 1 in the st position and 0's elsewhere. Then, since 



it is readily seen that 

E*E -e(e T e )e T - { ° iU ^ U 

-L^St-^UV c 8\ C t ^Uj^v ) Jjl _ u 



and hence that 



\2 



if s ^ t 
E st if s = t 



{Est, 

Let a = [ax, . . . , ai), ft = . . . , be multi-indices with 

(2. 1) positive integer entries, except for fa, which may also be zero 
and suppose further that 

(2.2) a 1 + --- + a e = a, fa + --- + fa = t, 
and let 

(2.3) m a)P (X, Y) = X ai Y h ■ ■ ■ X at Y Pt . 
Then, since s > £ and t > £ — 1 + fa, it follows that 

£ < s and £ < t + 1 - fa. 
The proof of the next two lemmas will rest heavily on the following observations: 
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if m, r and n are nonnegative integers such that m + r > 2 and n > 2, then 

[ E aa if a = b = c = d 
( 2 - 4 ) (£ cd ) m (£ a6 ) n (£ cd ) r = <{ o otherwise 

-Eaa if c = d = a 

(2.5) E aa E cd E aa = ^ o otherwise 

E aa if c = d = a 

(2.6) E ab E cc E ab = < q otherwise 

Remark 2.1. is a/50 useful to note that if the constraints (2.1) and (2.2) are in force 
and 

(2.7) X = XiE ab H and Y = XjE cd H , 

then 

(2.8) m Q>/3 (X, Y) = x^m^iEab, E cd ) + ■■■ . 

Lemma 2.2. Assume that the multi-indices a = (a l5 . . . , act) and (3 = . . . , /%) are 
subject to the constraints (2.1) and (2.2). Suppose further that s > 2, t > 2, and that 

(2.9) max{ai, . . . , a e ; . . . , /%} > 2. 
Then 

(2.10) majixiE^xjE*) = { * /a = J^T fl j • 

In other words, m a ^(xiE ab , XjE cd ) 7^ if and only if Xi and Xj are diagonal pairs in 
the same position. 

Proof. The proof is divided into cases. 

1. If £ = l,then 

m a)fi (xiE ab ,XjE cd ) = XiX^EahYi^E^Y 

and the asserted conclusion (2.10) follows from (2.4) with m = s and n = t, since 
s > 2 and t > 2, by assumption. 

2. If i > 1 and a r > 2 for some r e {1, . . . , fc}, then 

{x\E aa + ■ ■ ■ if b = a 
otherwise • 

But if 6 = a, then 

m^ixiE^xjEca) = x°x)(E aa )^(E cd )^E aa ) a2 • • • 
and (2.10) follows from (2.5). 
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3. If £ > 1 and ft > 2 for some r £ {1, . . . , then 

{x^-E cc if d = c 
otherwise • 

But if d = c, then 

m Q ^( Xl E ab , Xj E cd ) = X^Y^X^ ■■■ = z^$(£ a6 ) ai £ cc (£ o6 ) aa • • • 
and (2.10) follows from (2.6). □ 
Remark 2.3. Condition (2.9) is automatically met if either 

s > £, or t > £, or s = t = £ and ft = 0. 
It remains to consider the case 
(2.11) max{«i, . . .,a t ; ft, . . . ,ft} < 1. 

Lemma 2.4. If (2.1), (2.2) and (2.11) are in force and t > 2, then there are four 
possibilities: 

1. ft = 1: In this setting s = t, £ = s, m Qj/ g(X, Y) = (XYf and 

/ 77i \ s tfrp tp \t \ x\x\E aa ifc = bandd = a 

m a>p (x i E ab ,x j E ai )=x i x j {E ab E cd ) — | q otherwise ' 

2. ft = 0: In this setting s = t + 1, £ = s, ra a ^(X^ Y) = (XY fX and 

/ -w — i s4/rnrn . trn f xfx^Eab ifc = bandd = a 

m at p(x i E ab) x j E cd )=x i x j (E ab E cd ) E ab = j Q J otherwise ' 

3. ft = 1; /n zfe setting s = t, £ = s, m a ^(Y, X) = (FX)* and 



^(xjE^XiEah) - XiX^EdEabf - <j jj 1 ^ &a Z ^ C & Otherwise 



4. ft = 0: In this setting s = t + 1, £ = s, m a ,p(Y, X) = (YXfY and 

, th 7-i \ t 8/T-i T-i \t f XitfEbb ifc = bandd = a 
m a AxjE cd ,x t E ab ) = xlx^EcdEatf = j ^ /W/*? ' 

Proof. In view of (2.1), the constraint (2.1 1) implies that 

ai — ■ ■ ■ — on — 1, ft = • • • = ft-i = 1 and ft = 1 or ft = 0. 
Correspondingly 

f (XF) £ ifft = l 

m a ,p{X, Y) = 1 [xYf- l X if ft = 

and 

f (FX) £ ifft = l 

m 0hP (Y,X) = i (yxy^Y ifft = ■ 

The remaining conclusions are self-evident. □ 
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Definition 2.5. The two r letter monomials ex" 1 ■ • ■ xf r and fx^ 1 ■ ■ ■ x^ r with e ^ 

L\ L r <f Jl J r I 

and f are said to be > -equivalent if there exists a permutation a of the integers 
{1, . . . , r} such that j3j = d a ^) for j = 1, . . . r. This will be indicated by writing 

ex£---x? > /xf---xf . 

L\ L r " Jl Jr 

These two monomials are structurally equivalent (SE) if they are > -equivalent and 
e = f. This will be indicated by writing 

ex^---xf r SE /xj---xf. 
Thus, for example, if a, 6, c, d e M \ {0}, then the four two letter words 

24 r24 24 ^1^24 

CIX^Xqj OXqX^j CX^X^ HUQ (XX c^X q 

are [> -equivalent; they will be SE if and only if a — b — c — d. 
2.2. Enumerating one letter monomials in the k x k array. 

Lemma 2.6. If a family of polynomials pi, ■ ■ ■ ,Pk 2 in the 2k 2 commuting variables 
Xi, . . . , x 2 fc2 admits an nc representation p(X, Y), then for each positive integer n > 1 
exactly one of the following situations prevails: 

(1) There are no one letter monomials of degree n in any one of the given polyno- 
mials. 

(2) At least one of the given polynomials contains exactly one one letter monomial 
of degree n. 

(3) At least one of the given polynomials contains exactly two one letter monomials 
of degree n. 

Moreover, 

(2) holds <^=^ there exist exactly k polynomials each one of which contains exactly 
one one letter monomial ex™ of degree n (all with the same coefficient). 

(3) holds <^=^ there exist exactly k polynomials each one of which contains exactly 
two one letter monomials ex? + fx™" . 

L s ° Jt 

Further, if e ^ f and ef ^ then the letters Xi m , m = 1, . . . , k in the monomials 
exf l , . . . , ex™ are the diagonal entries of one of the matrices and the letters x Jm , m = 
I, ... ,k in the monomials fx^, . . . , /x™ fc are the diagonal entries of the other. 

Proof. Clearly (1), (2) and (3) are mutually exclusive possibilities that correspond to 

p(X, Y) = aX n + bY n + • • • 

with either (1) a = and b = for all n, (2) ab = and a ^ b 

Suppose first that at least one of the given k 2 polynomials contains exactly one term 
of the form ax™ with n > 1, a real coefficient a^O and Xj G X. Then 

X = XiE ss + ■ ■ ■ for some s e {1, . . . , k} 

and 

aX n = a(x™E ss + •••)• 
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Moreover, since the polynomial that contains the term ax n contains only one term of 
this form, it follows that 

p(X, Y) — aX n does not contain a term of the form cY n with c^O. 

Thus, as X has k diagonal entries, there will be exactly k polynomials each one of 
which contains a exactly one term of this form. This completes the proof of (a). The 
proof of (b) is similar to the proof of (a). 

Finally, a term of the form ax n + bx™ with i ^ j and ah ^ will be present in one 
of the polynomials if and only if either 

X = XiE ss + ■ ■ ■ and Y = XjE ss H 

for some choice of s £ {1, . . . , k}, or 

X = XjE ss H and Y = XiE ss H 

for some choice of s 6 {1, ... , k}. The rest of the proof goes through much as before. 

□ 

Lemma 2.7. Let pi,- ■ ■ , p^ be a family of polynomials with an nc representation 
p(X, Y). Suppose that some polynomial p in the family contains the terms ax n + bx™ 
with a b and either a ^ or b ^ 0. Then there exist k polynomials p^ , ■ ■ ■ , pi k 
containing the respective terms ax™ +&r™, ■ ■ ■ , ax™. +bx r - k . Moreover, the terms ax™ m 
with (1 < m < k) that have coefficient a are the diagonal variables of one matrix and 
the terms with bx™ m with (1 < m < k) that have coefficient b are the diagonal terms of 
the other matrix. 

Proof. This follows from Lemma 2.6. □ 

2.3. Enumerating one and two letter monomials in the k x k array. The symbol 

[ 1 if a y£ 

X(a) = < ifa = 

will be used in the next lemma. 

Lemma 2.8. Letpi, . . . , p^ be a family of polynomials in the 2k 2 commuting variables 
Xi, . . . , x 2 fc2 that admits an nc representation p(X, Y) of degree d > 1. Suppose further 
that s = t, t > 2, 

(2.12) p(X, Y) = e 1 (XY) t + e 2 {YX) 1 + e 3 X 2t + e 4 Y 2t + q{X, Y) 

where q is a polynomial that does not contain any scalar multiples of the first four 
monomials listed in (2.12), that the coefficients ei, . . . , are all distinct and that x u ^ 
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x v . Then the family ofk 2 polynomials will contain 



k 2 


- k 


2 letter monomials SE to eixtxt 

A it V 


if a ^ 


k 2 


- k 


2 letter monomials SE to e2X t u x t v 


ife 2 


k 2 


-k 


2 letter monomials SE to e^x^xl 


ife 3 ^ 


k 2 


- k 


2 letter monomials SE to e4X t u x t v 


ife 4 ± 


k 




2 letter monomials SE to <p(t, t)x^x* 


ifcp{t,t)^0 


k 




1 letter monomials SE to e$x^ 


ife 3 ± 


k 




1 letter monomials SE to e^x^ 


ife 4 



This list incorporates all the ways that one letter monomials of degree 2t and two letter 
monomials that are ^-equivalent to x^x* can appear in the given family of polynomi- 
als. Moreover, there is no cancellation: 

Each of the k polynomials that sit on the diagonal in the k x k array corresponding 
to p(X, Y) contains 

(x( e i) + X( e 2) + x( e 3) + x( e i))(k - 1) monomials > to a^x* 
made up of off-diagonal letters and 

(2.13) x{<p(ti t)) monomials > to x^x* 

made up of diagonal letters, as well as exactly one one letter monomial of degree 2t 
with coefficient e^, and exactly one one letter monomial of degree 2t with coefficient e$, 
both of which are diagonal entries. 

No off diagonal polynomial contains any two letter monomials > to x^x*. 

Proof. Two letter words x -x* with i ^ j and t > 2 may be generated in four different 
ways: 

(1) as entries in either m a ^{X,Y) or m a ^(Y,X) with a and (3 subject to (2.1) 
and (2.2) with s — t if Xi and Xj are in different matrices; 

(2) as entries in X 2t (resp., Y 2t ) if Xj and Xj are both in X (resp., Y). 

Suppose first that x« and Xj are in different matrices and that (2.9) is in force. Then 
Lemma 2.2 implies that the scalar multiples of the words x*x* will appear in at least one 
of the k 2 polynomials if and only if Xj and Xj are diagonal entries in the same position. 
In this instance, x\x^ will appear in a polynomial that sits in the same diagonal position 
as Xi and Xj with coefficient (f(t,t). 

Suppose next that Xj and Xj are in different matrices and that (2.1 1) is in force. Then, 
in view of Lemma 2.4, it remains only to consider the contributions from (17)' and 
{YXf: If 

X = XiE ab H and Y = XjE ba H , 

then 

(2.14) (XYf = x\x)E aa + ■■■ and (YX) 1 = x\x)E bb + ■■■ . 
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Since there are k 2 — k off-diagonal positions in a k x k matrix, there are k 2 — k choices 
of E ab with a ^ b. Moreover, since the entry x'x* appears in the aa position in (17)' 
and the bb position in (YX) 1 there will be no cancellation, even if e 2 = —e\. 

On the other hand contributions that come from diagonal entries of X and Y can 
interact with each other, i.e., if 

X = XiE aa H and Y = XjE aa H , 

then 

e^XYf + e 2 (YXY = ( ei + e 2 )x\x)E aa + ■■■ 

and 

e 1 {XY) t + e 2 (YXY + q{X, Y) = <p(t, t)x\x)E aa + ■■■ . 

Thus, if (p(t, t) 7^ 0, there will be k contributions, one for each choice of a G {1, ... , k}. 

The contributions from e^X 2t and e^Y 2t are enumerated in much the same way. 
Moreover, there is no cancellation, because the monomials with coefficient e 3 have all 
their letters in X and the monomials with coefficient e 4 have all their letters in Y . □ 

Remark 2.9. The list in Lemma 2.8 is written under the assumption that e\, e 2 , 
and are four distinct numbers. If, say, e\, e 2 and e 3 are three distinct numbers and 
e 4 = e3, then there will instead be 2(k 2 — k) two letter monomials SE to e^x^x^, k 2 ~k 
two letter terms SE to eix^x^, k 2 — k two letter terms SE to e^x^x^, k two letter terms 
SE to (p(t, tyx^xl and the one letter monomials would be as they are stated above. 

Lemma 2.10. Let pi, . . . ,Pk^ be a family of polynomials in the 2k 2 commuting vari- 
ables X\, . . . , x 2k 2 that admits an nc representation p(X, Y) of degree d > 1. Suppose 
further that t > 1, 

(2.15) p(X, Y) = hiXYfX + HYXfY + f 3 X 2t+1 + f,Y 2t+1 + q(X, Y), 

where q is a polynomial that does not contain any scalar multiples of the first four 
monomials listed in (2.15), the coefficients f\, . . . , are distinct and x u ^ x v . Then 
the family ofk 2 polynomials will contain exactly 



k 2 — k 2 letter monomials SE if fi ^ 

k 2 — k 2 letter monomials SE fax^x 1 ^ 1 if f 2 ^ 

k 2 — k 2 letter monomials SE j ^x^x^ i// 3 ^ 

k 2 — k 2 letter monomials SE f^x^x^ 1 if ^ 

k 2 letter monomials SE <p(t + 1, iftyif + 1, t) ^ 

k 2 letter monomials SE ip(t, t + if pit, t + 1) ^ 

k 1 letter monomials SE fzx^ +l if fz ^ 

k 1 letter monomials SE / 4 x^ t+1 if fi ^ 



This list incorporates all the ways that one letter monomials of degree 2t + 1 and two 
letter monomials > to can appear in the given family of polynomials. Moreover, 

there is no cancellation. 
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Each of the k 2 — k polynomials that are off the diagonal in the k x k array corre- 
sponding to p(X, Y) contains 

(x(fi) + XU2) + x(fa) + x(U)) two letter monomials > to x^xl 

made up of off-diagonal letters. Each of the k polynomials that are on the diagonal 
contains 

(2.16) x(y?(£ + 1, t)) + x(y?(t, * + 1)) two letter monomials > 

made up of diagonal letters, as well as a one letter monomial SE to fox^ 1 if fa 7^ 
and a one letter monomial SE to f±x^f +1 if f$ 7^ 0. The letters in these one letter 
monomials are diagonal entries. 

Proof. Let a = (oti, . . . , ai) and (/3i, . . . , fit) be multi-indices that meet conditions 
(2.1) and (2.2) and set s — t + 1. Two letter words with i 7^ j and t > 2 may 

be generated in four different ways: 

(1) by m ai p(X, Y) if Xi G X and Xj G Y; 

(2) by m a> p(Y, X) if X{ G Y and Xj G X; 

(3) by X 2t+1 if Xj G X and x 3 G Y; and 

(4) by Y 2t+l if Xj G F and Xj G Y. 

If (2.9) is in force, then the two letter words x\x l j with i ^ j can only come from 
diagonal pairs in the same position. The same terms with possibly different coefficients 
may appear from the diagonal entries in X and Y from polynomials of degree t + 1 in 
X and t in Y or degree i in X and degree £ + 1 in Y in q(X, Y) . The coefficients of the 
net contribution are ip(t + 1, £) and </?(£, t + 1), respectively, and there will be a total of 
kx(<p(t + 1, £) and t + 1) such pairs, one of each sort in each polynomial on 

the diagonal of the k x k array of p(X, Y). 
On the other hand, if (2.11) is in force, then 

«! = ••• = ocg = 1, /?! = ••• = = 1, = 

m Qi/3 (X,F) = (Xr)'X and m^(y, X) = 

Thus, if 

X = XiE ab H and Y — XjE cd H , 

then, in view of assertions 2 and 4 of Lemma 2.4, the coefficient of x\x^ in m a> p(X, Y) 
and m a; p(Y, X) will be nonzero if and only if c = b and d = a. Correspondingly, 

{XYfX = x l + l x)E ab + ■■■ and (YXfY = x\x] +1 E ba + ■■■ . 

Since there k 2 — k off-diagonal positions in a k x matrix, there are k 2 — k choices of 
i? a b with a ^ b. Moreover, the entry can not cancel the entry f2X t i x t J +1 even 

if a = b, since x^ and Xj are in different matrices. However, there can be contributions 
from monomials in q(X, Y) of degree t + 1 in X and t in Y or degree i in X and £ + 1 
in Y. 

Similarly, if a ^ b and X = XiE ab + XjE ba H (resp., F = XjE afc + £j-E fca H ), 

then 

X s+t = x\x)E ab + x*xj£; ba + • ■ ■ (resp., Y s+t = x s t x)E ab + x\x)E ba + ■■■). 
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The final assertion comes by counting the contributions discussed above. □ 

Remark 2.11. The list in Lemma 2.10 is written under the assumption that f\, f 2 , f'3 
and fi are four distinct numbers. If, say, f\, f 2 and f% are three distinct numbers and 
fi — f:i> then there will instead be 2{k 2 — k) two letter monomials SE to / 3 Xy +1 2r*, 
[k 2 — k) two letter terms SE to (k 2 — k) two letter terms SE to f 2 x t + 1 x t v , k 

two letter terms SE to <p(t + 1, and the one letter monomials would be as they 

are stated above. 

2.3.1. Enumeration of two letter monomials with one letter on the diagonal. 

Lemma 2.12. Let pi, ■ ■ ■ ,Pk 2 be a family of polynomials in the 2k 2 commuting vari- 
ables xi, . . . , x 2k 2 that admits an nc representation p(X, Y) of degree d > 1. Suppose 
further that t > 2, 

(2.17) p{X, Y) = diX 1 + d 2 X t - 1 Y + d^YX 1 - 1 + q(X, Y), 

where q(X, Y) does not contain any scalar multiples of the first three monomials in 

(2.17) and the coefficients di, . . . , d% are distinct. 

If x u is a diagonal element of X, then the family of k 2 polynomials will contain 
exactly 

2k — 2 2 letter monomials d\x l ~ x x v with x v 7^ x u if d\ 7^ 

k — 1 2 letter mnmls d 2 x t ~ l x v with x v an off-diagonal entry ofY if d 2 7^ 

k — 1 2 letter mnmls d^x l ~ l x v with x v an off -diagonal entry ofYifd^^O 

1 2 letter mnml (p(t — 1, with x v a diagonal entry ofY if(p(t — 1, 1) 

This list incorporates all the ways that two letter monomials of degree t with x u 
of degree t — 1 can appear in the given family of polynomials. Moreover, no two 
monomials in this list are the same. 

If x u is in the aa position of X and p a b denotes the polynomial in the ab position in 
the k x k array corresponding to p(X, Y), then 

(2.18) pa^xx,. . . ,x k 2) = dix l ~ l x v + d 2 x t ~ 1 x w ^ for a ^ b, 

where x v (resp., x w ) is in the ab position in X (resp., Y) and there are no other two 
letter monomials of degree t in p ab with x f ~ x as a factor. Similarly, 

(2.19) Pba(xi, ■ ■ • ,^fc 2 ) = diX^Xm + d 2 x t ~ l x n H for a ^ b, 

where x m (resp., x n ) is in the ba position in X (resp., Y) and there are no other two 
letter monomials of degree t in p ba with x l ~ x as a factor. 

Proof. It is readily checked with the aid of the calculations in §2.1 that if 

(1) x u is in the aa position of X, x v E X and x v 7^ x u , then x l ~ x x v is either in the 
ab position of X 1 or the ba position of X 1 for some b 7^ a; 

(2) x u is in the aa position of X and x v E Y, then x t ~ 1 x v is either in the ab position 
of X l ~ l Y or the ba position of YX l ~ x for some b 7^ a. 

The rest of the proof is straight forward counting and is left to the reader. □ 
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Remark 2.13. Let 

(2.20) 

p{X, Y) = diX* + d 2 X t - 1 Y + dsYX 1 - 1 + rf 4 y* + q(X, Y) for some integer t > 2, 

where q(X,Y) does not contain any scalar multiples of the first four monomials in 

(2.20) and assume that dx, . . . , are subject to the constraints 

(2.21) di 7^ 0, di / <i 4 and d\ ^ d 2 or d\ ^ d 3 . 
Then there will be k terms 

d\x^ , . . . , d\x ik 

in the family ofk 2 polynomials and the corresponding letters x\ v . . . , x\ may be iden- 
tified as the diagonal elements of say X. The assumption d\ ^ insures that they can 
be chosen unambiguously. Thus, if x u is one of these diagonal elements and it is in the 
aa position of X, then 

Paa = dxx* u + <f(t - 1, l)x l ~ l x h H , 

p ab = dix t ~ 1 x v + d 2 x t ~ 1 Xf H , 

Pba — d\X l u 1 X Z + d^X l u 1 X g + • • • , 

where Xh is in the aa position ofY, x v is in the ah position ofX, x / is in the ah position 
ofY, x z is in the ha position of X and x g is in the ha position ofY. 

2.4. Determining Diagonal Elements. Lemmas 2.8 and 2.10 serve to enumerate the 
diagonal entries in X and Y when the given set of k 2 polynomials contain one letter 
monomials. But if say e% = in Lemma 2.8 and fs = in Lemma 2.10, then it is 
not immediately obvious which entries belong to X and which entries belong to Y. 

Definition 2.14. A pair of variables X; L and Xj will be called a partitioned (resp., 
dyslexic) diagonal pair if both occur in the aa position for some choice of a £ 
{1, . . . , k} and we know (resp., do not know) which variable occurs in X and which 
occurs in Y. 

Lemma 2.15. Suppose Pi, ■ ■ ■ ,p 2 k 2 is a family of polynomials in the 2k 2 commut- 
ing variables x\, . . . , x 2 ^ that admits an nc representation p(X, Y) of degree d with 
d > 1 such that at least one of the given polynomials contains a term of the form 
exfm(xi, ...,Xi, x 2 fc2), where eeK \ {0} and m(xi, ...,Xi, x 2 fc 2 ) is a monomial 
of degree t > 1 that does not contain any Xi terms and s > t + 2. Then Xi lies on the 
diagonal of either X or Y . 

Proof. For the sake of definiteness, assume Xi £ X. Then, since s > t + 2, every 
permutation of the symbols X s Y l must contain at least two adjacent X's. But if 

X = XiE cd + ■ ■ ■ with d, 

then X 2 = 0. Thus, the given family of polynomials will only contain terms of the 
form exfm(xi, Xi, x^) if c = d, i.e., if X{ is a diagonal element of X. □ 

Lemma 2.16. Let pi, . . . ,Pk^ be a family of polynomials in the 2k 2 commuting vari- 
ables xi, . . . , x^ that admits an nc representation p(X, Y) of degree d > 1. Suppose 
further that at least one of the polynomials contains at least one term of the form ex^x^ 
with s > t + 2, and t > 2. Then 
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(1) Xi and Xj are a dyslexic diagonal pair. 

(2) There exist exactly k dyslexic diagonal pairs {x^, x^}, . . . , {xi k , Xj k } and k 
polynomials p^, . . . ,pi k such thatp in = ex^x^-l forn = 1, . . . , k. More- 
over, pi n occupies the same diagonal position as x in and Xj n . 

Proof. If Xi E X, then (since s > t + 2 and t > 2) Lemma 2.15 guarantees that Xi 
lies on the diagonal of X and that Xj lies on the diagonal of the matrix that contains 
Xj. Since Xi and Xj must be in the same diagonal position, this forces Xj to belong 
to Y. Similarly, if Xi E Y, then Xj must be in X and both variables must be in the 
same diagonal position. Thus (1) holds; (2) follows automatically from (1), since all 
the dyslexic diagonal pairs are subject to the same constraints. □ 

Lemma 2.17. Let pi, ■ ■ ■ ,Pk 2 be a family of polynomials in the 2k 2 commuting vari- 
ables xi, . . . , x 2 fc2 that admits an nc representation p(X, Y) of degree d > 1. Ift > 2 
and s = t + lors = t and it is also assumed that 

(1) k > 3 (so that k 2 — k > k) and 

(2) there exist exactly k monomials exf^x^, ■ ■ ■ , exf ' Xj in the given family of 
polynomials that are structurally equivalent to ex s u x l v with x u ^ x v and e ^ 0, 

then x in and xj n are a dyslexic diagonal pair for each 1 < n < k. Moreover, if 
Pe„ = ex l n xt j n + • • • » then pe n occupies the same diagonal position as Xi n and Xj n . 

Proof. Since there are only k terms in this list and k 2 — k > k when k > 3, the 
conclusion follows from Lemma 2.8 if s = t and from Lemma 2.10 if s = t + 1. In 
the first case e = cp(t, t); in the second case either e = cp(t + 1, t) or e = cp(t, t + 1) 
and the two numbers cp(t + 1, t) and cp(t, t + 1) are different. □ 

Lemma 2.18. Suppose that the term ax\x^, a / 0, occurs in some polynomial p in the 
given family. If the pair Xi and Xj is a dyslexic diagonal pair, then the polynomial is a 
diagonal polynomial that occurs in the same diagonal position as the dyslexic pair. 

Proof. If, say, X = XiE aa + • • • , and Y = XjE aa + ■ ■ ■ and the multi-indices a and (5 
are subject to (2.1) and (2.2), then m Q , i/3 (X, Y) = x\x l jE aa + • • ■ . □ 

Lemma 2.19. Let pi, ■ ■ ■ ,P2k 2 be a family of polynomials in the Ik 2 commuting vari- 
ables X\, . . . , X2k 2 that admits an nc representation p(X, Y) of degree d with d > 1 
such that exactly k SE terms of the form ax^x^, ■ ■ ■ ,ax| fe a;* fe appear in the family 
with s > t, t > 2 and a / 0. If '{x^, x^}, ■ • ■ , {x ik ,Xj k } are dyslexic diagonal pairs, 
then the variables x^, • • • , Xi k of degree s in the terms ax^x^, ■ ■ ■ , ax^x^ are the 
diagonal elements of one matrix and and the variables Xj 1 , ■ ■ ■ , Xj k of degree t are the 
diagonal elements of the other matrix. 

Proof. If the lemma is false, then without loss of generality, we may suppose that 
E X and x i2 E Y. By Lemmas 2.16 and 2.10, p(s,t) = <p(t,s) = a ^ 0. Thus, 
as s > t, each diagonal pair X; , Xj will occur in the monomials ax? x\ and ax] X s , . 
Therefore, there will be 2k terms structurally equivalent to ax^x^, which contradicts 
one of the given assumptions. □ 
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3. Partitioning Algorithms for families containing single letter 

MONOMIALS ax?, n>2 

In the previous section we developed a number of methods to determine the diago- 
nal variables for a family V of polynomials p 1 , ■ ■ ■ ,p k 2 with an nc representation. The 
next step is to utilize this information to determine which of the commutative variables 
• , x 2 k 2 are entries in X and which are entries in Y. In this section we develop 
two algorithms that we will refer to as partitioning algorithms, since they partition 
the commutative variables between the two matrices. First, however, we review some 
preliminary calculations that will be essential in the development of the first partition- 
ing algorithm. 

The following assumptions will be in force for the rest of this section: 

(Al) pi, . . . ,p k 2 is a family of polynomials in 2k 2 commuting variables x\, . . . , x 2 k 2 
that admits an nc representation p(X, Y) of degree d with d > 1. 

(A2) ip(n, 0) 7^ for some n > 2. 

(A3) p(0, 0) = 0. (This involves no real loss of generality, because constant terms 
can be reinserted at the end.) 

3.1. Preliminary calculations for the partitioning algorithm. First observe that 

(aI k + /3E st ) n = a n I k + na n - l (3E st lis^t. 

Then without loss of generality, assume that the variables {x ii: ■ ■ ■ , x ik } are the com- 
mutative variables that occur on the diagonal of X. 
Next, choose a variable Xj such that 

Xj ^ {2;^ , . . . , Xi k }, 

set 

Xjj = • • • = x ik = a and Xj = (3 

and set all the other variables equal to zero. 

Then there are three mutually exclusive possibilities: 

1. Xj E X. In this case 

X = al k + f3E st with s ^ t and Y = 0. 

2. Xj G Y but is not on the diagonal of Y. In this case 

X = al k and Y = /3E st with s^t. 

3. Xj G Y and is on the diagonal of Y. In this case 

X = al k and Y = /3E SS for some integer s G {1, k}. 



22 H. DYM, B. HELTON, AND C. MEIER 

All three of these cases fit into the common framework of choosing X = A e 
and Y = B e R kxk with AB = BA. Thus, if 

d 

p(X,Y) = Y,PdX,Y), 

i=0 

where P[{\{X, Y) denotes the terms in p(X, Y) of degree i, then the condition AB 
BA insures that 

i i 

(3.1) Pli] (A, B) = Y, * j.j-V '!>> ! = + ]T <; ,,.1' 'II' = p K (A, B), 

j=0 j=l 

where 

Cij = (p(i,j) (for short) for i, j — 0, . . . , d with Qj = for % + j > d 

and hence, 
(3.2) 

d i d d i 

P (a, B) = J2J2 '•< j-j a ' = E c * Al + E E <■ j-j a ' j & = p(a b), 

i=l j=0 i=l i=l j=l 

since the assumption p(0, 0) = forces c 00 = 0. 

In Case 1, A = al k + /3E st with s and B = 0. Therefore, 

(3.3) P\i](A,B) = Cioipth + ia^PEst) fori > 1 and 

d d 

p(A,B) = E^o^ + E ^' 1 ^ 
i=i i=i 

(3.4) = e c *« l 4 + E + ^v^- 

j=l i=0 

In Case 2, A = al k and £? = (5E st with s^i. Therefore, 

(3.5) P[i\(A B) = c i0 a l I k + Ci^ hl a l ~ l f3E st for % > 1 and 

(3.6) p(A, 5) = E c *o«^ + E 

i=i i=i 

In Case 3, A = and £? = /3-E ss . Therefore, 

« 

(3.7) B) = c i0 a l I k + E c^jja 1 ' 3 (3 J E ss for i > 1 and 

i=i 

(3.8) p(al k , (3E SS ) = <*)<*% + E E E S8 . 
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We remark that the formula p(alk + /3E st , 0) in Case 1, can also be expressed in 
terms of the polynomial 

d 

<P(<X) = ^CiQO 1 
i=l 

as 

p(al k + (3E st , 0) = <p{a)I k + tp'(a)pE st . 
The equality q = a ^ guarantees that (p(ct) ^ and ip'(a) ^ 0. 

3.2. Partitioning Algorithm I. Now we develop the fundamental ideas for an algo- 
rithm that will partition the commutative variables between X and Y. 

Set the diagonal entries of X equal to a, one of the other 2k 2 — k variables xi equal 
to (5 and the remaining 2k 2 — k — 1 variables equal to zero. Then: 

Xi is an off-diagonal entry of X -<=>- this substitution produces 

d 



k polynomials equal to ^ c^ol\ 



i=i 

d-l 



1 polynomial equal to ^ Q+i,o(^ + and 
k 2 — k — 1 polynomials equal to 0; 



Xi is an off-diagonal entry of Y -<=>- this substitution produces 

d 

k polynomials equal to ^ c i0 a l , 



i=i 

d-l 



1 polynomial equal to ^ Cna l fi and 



i=0 



k 2 — k — 1 polynomials equal to 0; 



Xi is a diagonal entry of Y -<=>- this substitution produces 

d 

k — 1 polynomials equal to ^ c i0 a\ 

i=\ 

d d i 

1 polynomial equal to ^ c i0 a l + ^ ^ Ci^jja l ~ j (3 j and 

1=1 j=l jr' = l 

k 2 — k polynomials equal to 0. 
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The first two cases will be indistinguishable if and only if 

d-l d-l 

Ci + i >0 (i + l)a\0 = Ciiai P 

i=0 i=0 

for every choice of a and f3, i.e., if and only if 

Ci+i,o(i + 1) = Cji for i = 0, . . . ,d - 1. 

3.3. Partitioning with homogeneous components: Algorithm DiagParl. It is often 
advantageous to focus on the homogeneous components of the given set of polynomials 
Pi, . . . ,p k 2, i.e., on the sub-polynomials of p 1: . . . ,p k 2 of specified degree. It is readily 
seen that if (X, Y) denotes the terms in the nc polynomial of degree n, then the 
k 2 commuting polynomials in the k x k array corresponding to p[ n ] (X, Y) are the 
k 2 polynomials qi, . . . , q k 2, where q^xi, . . . ,x 2k 2 ) is the sum of the monomials in 
Pi(xi,.. . , x 2k 2) of degree n. 

The three possibilities considered earlier applied to polynomials of degree n lead to 
simpler criteria: 

1. If A = al k + f3E st with s ^ t and B = 0, then 

p [n ](A,B) = c n0 {a n I k + na^pEst). 

2. If A = al k and B = (3E st with (3 ^ and s ^ t, then 

p [n] (A,B) = c n0 a n I k + Cn-^a^pEst. 

3. If A = al k and B = (3E SS , then, 

n 

P[n](A, B) = c n0 a n I k ■ /•/"" J ;,/ '-'- 

3=1 

Thus, if the diagonal entries of X are set equal to a, one of the other 2k 2 — k variables 
Xi is set equal to (5 and the remaining 2k 2 — k — 1 variables are set equal to zero, and 
if (5 7^ a and ip(n, 0) ^ 0, then: 

Xi is an off-diagonal entry of X <^=^ this substitution produces 

k polynomials equal to c n0 a n , 

1 polynomial equal to c n0 na n ~ 1 /3 and 

k 2 — k — 1 polynomials equal to 0; 



Xi is an off-diagonal entry of Y <^=^> this substitution produces 

k polynomials equal to c n oa n 

1 polynomial equal to c n _i i ia: ri_1 j3 and 

k 2 — k — 1 polynomials equal to 0; 
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Xi is a diagonal entry of Y -<=>• this substitution produces 
k — 1 polynomials equal to c n oa n 

n 

1 polynomial equal to c n0 o; n + 2_, c n -jja n ~ j (3 j and 

k 2 — k polynomials equal to 0. 

(1) The first two possibilities will be distinguishable if and only if 

nip(n, 0) 7^ ip{n — 1,1). 

(2) The third will be distinguishable from the first if Lp(n, 0) 7^ by counting the 
number of polynomials equal to zero, k + 1 vs k 

(3) The third will be distinguishable from the second if and at least two of the 
coefficients tp{n — j,j) 7^ j = 0,1, ... ,n are nonzero. This is done by 
comparing number of terms of polynomials. 

We will refer to the process developed in the previous discussion as Algorithm 
DiagParl. In the following theorem we summarize the conditions under which this 
algorithm will partition the commutative variables. The reader should keep in mind 
that if there are k single letter monomials (as opposed to 2k or 0) of degree n in a 
homogeneous family V of polynomials of degree n, then we always assume that the 
associated variables lie on the diagonal of X and hence that ip(n, 0) ^ 0. 

Theorem 3.1 (DiagParl Algorithm). Letpi, ■ • ■ , p k 2 be a family of homogeneous poly- 
nomials with an nc representation p(X, Y) of degree n with n > 2. Then Algorithm 
DiagParl successfully partitions the variables x\, ■ ■ ■ , x 2 k 2 between X and Y if and 
only if 

(3.9) (p(n, 0)^0, ip{n,0) ^ ip{0,n), and n<p(n, 0) <p(n — 1,1). 

Proof. If cp(n, 0) 7^ and (3.9) holds, then the above discussion implies that DiagParl 
will successfully partition the variables between X and Y. 

Conversely, suppose that DiagParl successfully partitions the variables between X 
and Y. This implies that the algorithm can first determine the set of diagonal variables 
of X by analyzing single letter monomials of the form ax™ in the polynomials in V . 
By Lemma 2.6, we see that this is possible only if Lp(n, 0) 7^ and (p(n, 0) 7^ ^(0, n). 
Finally, the above discussion implies that DiagParl partitions the off-diagonal elements 
between X and Y if and only if tp(n, 0) 7^ tp(n — 1,1). □ 

3.4. Partitioning: Algorithm DiagPar2. In this subsection we develop another par- 
titioning algorithm that is closely related to theDiagParl Algorithm and is based on the 
following observation: Let Pi, • • • ,Pk* be a family of polynomials with an nc repre- 
sentation p(X, Y) and suppose that for some t > 2, (p(t, 0) 7^ and (p(t, 0) 7^ y?(0, t). 
Then by Lemma 2.6, there exists exactly k polynomials Pi x , - ■ ■ ,Pi k , each of which 
contains exactly one one letter monomial of degree t with coefficient cp(t, 0). 
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The next step rests on Lemma 2.12. But for ease of understanding, let xij denote 
the ij entry in X, y^ the ij entry in Y and let denote the ij entry in an array of 
commutative polynomials that admits an nc representation p(X, Y). Then 

Pu = <p(t,0)x t H + (p(0,t)y t u + -- - , 

whereas for i ^ j, 

Pij =<p(t, 0)(x*r 1 Xij + .r, r r l jj 1 ) + p(0, t){y\r l yij + y^f) 
+ <p(Y; t - 1, 1 )y ij x% : + <p(t - 1, 1; Y)x t rr l y i] + ■■■ 

and 

Pji =(p(t, 0)(x ii x*r 1 + xffxji) + <p(0, ^{yjiyli 1 + yjj 1 ^) 
+ if(Y; t - 1, \) yjl x\- x + <p(t - 1, 1; YU'j/uy, + • • ■ 

Thus, if the entry Xjj and £ are known, then there will be exactly two two letter monomi- 
als in p^, i ^ j with x l u l as a factor, namely, <p(t, Q)xj i x t i ^ 1 and cp(t — 1, 1; Y)x t ^ l y i j, 
but only the first of these will have the correct coefficient if (p(t, 0) 7^ cp(t — 1, 1; Y). 
Thus, under this condition, it is possible to isolate all the entries in X by repeating the 
argument for i = 1, . . . , k. 

Similar considerations based on inspection of the polynomials p^ will also yield all 
the entries in X if <p(t, 0) ^ <p(Y; t - 1, 1). 

We will refer to the above procedure for partitioning the commutative variables as 
Algorithm DiagPar2. The conditions under which this algorithm works are summa- 
rized in the following theorem. 

Theorem 3.2 (DiagPar2 Algorithm). Let Pl , ■ ■ ■ ,p k 2 be a family of polynomials with 
an nc representation p(X, Y). Then Algorithm DiagPar2 will partition the variables 
x ii ' ' ' ? x 2k 2 between X and Y if and only if there exists an n > 2 such that <p{n, 0) 7^ 
0, <p(n, 0) 7^ (p(0, n), and 

(3.10) (p(n,0) ^cp(n- 1,1; Y) or cp(n, 0) ^ (p(Y;n - 1, 1). 

Proof. If the above conditions hold then the above discussion implies that Algorithm 
DiagPar2 will successfully partition the variables between X and Y. Conversely, if 
DiagPar2 partitions the variables between X and Y, then it must first determine the 
diagonal elements. Lemma 2.6 implies that this is only possible if (p(n, 0) 7^ and 
<p(n, 0) 7^ cp(0, n). The above discussion implies that DiagPar2 will successfully par- 
tition the off-diagonal entries only if the above conditions hold. □ 

3.5. Summary of Partitioning Algorithms. The main conclusion of this section is 
that if the given system V of k 2 commutative polynomials admits an nc representation 
p{X, Y) in the set NC^xi) of nc polynomials p{X, Y) of degree d > 1 for which there 
exists an integer n > 2 such that 

{(p(n, 0) 7^ 0, ip(n, 0) 7^ ip(0, n), and either 
n<p(n, 0) ^ tp{n- 1,1) or 
(p{n,0)^(p{n-l,l;Y) or <p(n, 0) ^ (pQf; n - 1, 1), 

then either Algorithm DiagParl or Algorithm DiagPar2 will successfully partition the 
variables between X and Y, 
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4. Positioning algorithms for families of polynomials containing 

one letter monomials 

In this section we shall present algorithms for positioning the variables in X, given 
that the diagonal entries of X are known and that the remaining 2k 2 — k commutative 
variables are partitioned between X and Y. The assumptions (Al), (A3) that are listed 
at the beginning of §3 and a weaker form of (A2): 

(A2') \(p(n, 0) | + \ip(0, n)\ > for some n > 2, 

will be in force for the rest of this section. 



4.1. Positioning the variables within X: Algorithm ParPosX. If we can determine 
the diagonal variables of the matrix X and implement Algorithm DiagParl or DiagPar2 
in §3, then we may assume that x^, . . . ,x ik2 £ X, and the remaining k 2 variables 
belong to Y. To ease the notation, assume that Xi, . . . , x k 2 £ X and that Xi is in the ii 
position for i — 1, . . . , k and let 

Ri denote the remaining k — 1 entries in the ith row of X, 

Ci denote the remaining k — 1 entries in the ith column of X, and 

Li = Ri U Ci. 
Then, since 

(4.1) X n = (x^Eu + ■■■ )( Xj E st + ■■■) = (x^XjEuEst + ■■■) 
and 

(4.2) X n = ( Xj E st + ■■■ + •••) = (x^xjE^Eu + ■■■) 

for 1 < i < k < j, it is readily seen that the term ax^'^Xj appears in one of the entries 
of aX n if and only if either s = i or t = i, i.e., if and only if Xj £ Ri U Cj. Moreover, 
since Ri fl Ci = there will be 2k — 2 such terms in X n . 
Let 

Li fl L r = {xj, x e } for some integer r £ {1, . . . , k} \ {i}. 

Then one of these two variables will be in the ir position of X, while the other is in the 
ri position, and it is impossible to decide which is where without extra information. 
Let us assume for the sake of definiteness that Xj is in the ir position of X, then xi will 
be in the ri position. Moreover, since 

X n = ( Xj E ir + ■ ■ -){x n r - l E rr + ■■■) = (xjx^Eir + •••), 

the term XjX™~ 1 also belongs to the same polynomial. Thus, 

"<?n ■ ■ ■ Qik 

X n = : : 

.Qkl ■ ■ ■ Qkk. 

where the q st are either homogeneous polynomials of degree n in the variables x\ , . . . , x k 2 
or zero. In particular, if k > 3, n > 2 and i ^ r, then 

Qir ip^i %j%r ^i %t%m ~t~ ' ' ' ) 
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with x t and x m off the diagonal and x t ^ x m . This automatically insures that x t and 
x m differ from Xi and Xj, i.e., 

{x t , x m } n {x h xj} = 0, 

and that there exists an integer j G {1, . . . , k} \ {i, r} such that 

either x t G Ri D Cj and x m G -Rj D C r or vice versa. 

Therefore, since Ri DRj = for j ^ i, only one of the two variables x t , x m is listed in 
the set Ri U Cj and hence these two variables can be positioned unambiguously. This 
procedure adds one more variable to each of Ri and C r . The remaining entries in Ri 
and C r are obtained by repeating this procedure k — 3 more times by running through 
all the other triples of the form x™~ 2 x u x v in q ir . 

After Ri and C r are filled in, the procedure is repeated with some other diagonal 
element as a starting point, and then repeated again and again until all the k 2 variables 
are positioned in X. 

We will refer to the procedure for positioning the variables outlined in the previous 
discussion as Algorithm ParPosX. The conditions under which this algorithm will 
position the commutative variables are summarized in the following theorem. 

Theorem 4.1. [Algorithm ParPosX] Let p%, - ■ ■ ,Pk 2 be a family of polynomials with 
an nc representation p(X, Y). Then Algorithm ParPosX will position the variables in 
X if and only ifp G iVC^.n). 

Proof. DiagParl and DiagPar2 will partition the commutative variables between X 
and Y if and only if p G NCt s .ii)- The above algorithms will position the variables in 
X if and only if the commutative variables are partitioned between X and Y. □ 

4.2. Algorithm for positioning the polynomials given positioning in X. If X and 

Y are general matrices containing the variables x%, ■ ■ ■ ,x 2 k 2 , an nc representation 
p(X, Y) produces a family of polynomials Pi, • • • ,Pk 2 that is arranged in a k x k 
matrix. The purpose of this section is to investigate when it is possible to determine 
the position of the polynomials pi, • • • , Pk 2 in the resulting kx k matrix. We start with 
a lemma that will be useful in developing a procedure to accomplish this task. 

Lemma 4.2. Let pi, ■ ■ ■ ,p k 2 be a homogeneous family of polynomials with an nc rep- 
resentation p(X, Y) of degree d > 1. Suppose that the variables x%,--- , x 2 k 2 have 
been partitioned between X and Y and positioned in X. Then, if Xi is in the aa posi- 
tion of X, Xj is in the ah position of X, xi is in the bb position of X and 

ip(n, 0) 7^ for some n>2, 

then there exists exactly one polynomial that contains the monomials <p(n, 0)x™~ 1 Xj 
and (p(n, 0)xjX™~ 1 . Moreover, this polynomial is in the ab position in the k x k array. 

Proof. This is an easy consequence of the discussion in Subsection 4.1. □ 

Lemma 4.2 suggests a very simple algorithm for positioning the polynomials in our 
family. If the commutative variables are partitioned and positioned in the matrix X 
and if (p(n, 0) ^ for some n > 2, we simply run through the terms cp(n, 0)x™~ 1 Xj 
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for each X{ in the ii position in X and Xj in the ij position in X. The polynomial that 
contains this monomial will be in the ij position in the array of polynomials. We will 
refer to this procedure as Algorithm PosPol. The next proposition summarizes the 
conditions under which it is applicable. 

Proposition 4.3 (Algorithm PosPol). Letpi, • • • , p k 2 be a homogeneous family of poly- 
nomials with an nc representation p(X, Y) of degree d > 1. Then Algorithm PosPol 
will successfully position the polynomials p±, • • • , p k z if and only ifp e NC^ny 

Proof. Algorithm PosPol will position the polynomials if and only if the 2k 2 variables 
are partitioned, the variables in X are positioned and cp(n, 0) ^ 0. Algorithm ParPosX 
will position the variables in X if and only ifpE NCr^.n)- □ 

4.3. Positioning Y given the position of the polynomials: Algorithm PosY. Given 
a family V of polynomials p±, • • • ,p k 2 in the variables Xi, • • • ,x 2 k 2 with an nc rep- 
resentation p(X, Y), we have to this point developed algorithms to partition the vari- 
ables between X and Y, position the variables in X and position the polynomials 
Pi, ■ ■ ■ ,p k 2. The next step is to develop an algorithm that positions the commutative 
variables in Y. 

Suppose that the commutative variables x 1 , - ■ ■ , x 2k 2 nav e been partitioned between 
X and Y and positioned in X. Furthermore, suppose that the polynomials in the 
family V have been positioned and that cp(s, t) ^ for some s > and t > 1. Let pij 
denote the polynomial in the ij position in p(X, Y). Let x\, ■ • ■ , x^ be the variables 
contained in X and let X\, ■ ■ ■ ,x k be the diagonal variables of X. Set x\ = ■ ■ ■ = 
Xk = 1 and Xk+i = ■ ■ ■ = x^ = 0. On the nc level this is equivalent to setting 
X = I k . Then consider the polynomials %j = Pij(l, ■ ■ ■ , 1, 0, • • • ,0, x k 2 + x, ■ ■ ■ , x 2 k 2 ) 
in k 2 commuting variables. If V has an nc representation p(X, Y), this collection of 
polynomials will have an nc representation p(I, Y) = p{Y) that contains the monomial 
cp(s, 7^ 0. Therefore, each diagonal polynomial pa will contain a monomial of the 
form ip(s, where ji > k 2 + 1. This implies that Xj t must be the ii entry of Y. By 
repeating this argument for each i we can position the other diagonal elements of Y. 
To position the remaining variables, observe that a monomial of the form (p(s, t^x^Xu 
will appear in each %j where u > k 2 + 1. Given that Xj t lies in the ii position of Y 
and pij is in the ij position in p(Y), it follows that x u is in the ij position in Y. By 
repeating this argument for each ij we can position the other non diagonal entries of 
Y. 

We will refer to the process described above as Algorithm PosY. We summarize the 
conditions under which it will successfully position the variables in Y in the following 
proposition. 

Proposition 4.4. [Algorithm PosY] Let pi, • • • ,Pfc 2 be a family of polynomials with 
an nc representation p(X, Y). Then Algorithm PosY will successfully position the 
commutative variables in Y if and only if p G NC@ .in and there exists a pair of 
integers s > and t > 1 such that (p(s,t) ^ 0. 



Proof. This follows from the preceding discussion and the fact that PosPol will posi- 
tion the polynomials if and only if p 6 NC^.n). □ 
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4.4. Uniqueness results for one-letter algorithms. In this section we investigate the 
possible variations in the matrices X and Y that are obtained by Alg.l and Alg.2. We 
shall assume that the algorithms are applied to the homogeneous components of degree 
n in the family of polynomials V, where n > 2 and is also the lowest degree of single 
letter monomials in V. 

If P denotes an array of the given k 2 polynomials (as on the right hand side of (1.3)) 
that admits the nc representation 

(4.3) P = p(X, Y) = ^ c a> i3m a _p(X, Y), where \a\ + \/3\ = n in the sum 

and II is a k x k permutation matrix, then 

(4.4) n T Pn = c a ,pm a .^U T XU, U T YU). 

Moreover, since 
it is readily seen that 

(4.5) P r = ^c^m /3V ,(r T ,X T ), 
where 

P = (ft, ...,p r )=>p'= (p rt ...,Pt) and a = (a u . . . , a r ) ^> a' = (a r , . . . , an). 

Below we shall show that, aside from a possible interchange of X and Y, formulas 
(4.4) and (4.5) account for the only possible variation in the matrices X and Y that are 
generated by our algorithms. They correspond to the fact that the ParPosX algorithm 
allows for the diagonal variables to be placed in any order along the diagonal of X 
and the ambiguity in the next step of that algorithm, in which two variables x u and 
x v are arbitrarily assigned to be either the st entry or the ts entry of X (for some non 
ambiguous choice of s and t with s ^ t). This is in fact the only freedom that one has 
in positioning the variables within X, i.e., once these k + 2 variables are allocated, the 
positions of the remaining variables in X are fully determined. This is substantiated 
by the next theorem. 

Recall that the first matrix constructed by either Alg.l or Alg.2 is always designated 
X(or X). 

Theorem 4.5. Suppose that V is a homogeneous family of polynomials Pi, • • • , Pk 2 of 
degree n in the commutative variables x%, • • • , X2k 2 that admits an nc representation. 
Then if X, Y and X, Y are generated by two different applications of Alg.l or Alg.2, 
they are permutation equivalent (as defined in (1.13)). 

Proof. Suppose first that the given family V of k 2 polynomials contains exactly k 
one letter monomials (all with the same coefficient a G E \ {0}) and 

let X and X denote the matrices that are determined by successive applications of 
either Alg.l or Alg.2. Then the set of diagonal entries of X (without regard to to their 
positions on the diagonal) will coincide with the diagonal entries of X. Moreover, 
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since D2Pal and D2Pa2 depend only upon the diagonal variables and not upon how 
they are positioned along the diagonal, the set of k 2 variables in X coincides with the 
set of k 2 variables in X. 
Let 

L(x ia ) = {x u G X \ {x is } : ax" _1 x u appears in one of the polynomials in V}. 

There are 2k — 2 distinct variables x u in L(x is ). Moreover, the selection of these terms 
is based totally on V and not upon the position of x is in X. Therefore, if s ^ t, then 
the two variables in the intersection 

L(xi s ) H L(x it ) = {x u , x v }, 

are also independent of the position of x is and x it . 

If x is is the ss entry of X and x it is the tt entry of X, then either x u is the st entry 
and x v is the ts entry, or vice versa. It is impossible to decide. However, once one 
of these two possibilities is chosen, then the position of all the remaining k 2 — k — 2 
variables in X are determined by the algorithm. 

The algorithm ParPosX accomplishes this by inspecting terms of the form x™~ 2 x t x m 
in V, and so once one variable is placed, the algorithm is able to partition the set L(x is ) 
into the sets 

R(x is ) = {x u G X : x u is in the same row as x is } 

and 

C(x is ) = {x u G X : x u is in the same column as x is } 

for each integer s, 1 < s < k, solely by analyzing monomials contained in V. Con- 
sequently, these sets do not depend on the position of the diagonal entry Xi g . Since 
X has the same diagonal elements as X, there exists a permutation n of the integers 
{].,...,&} such that if x ss denotes the ss entry in the matrix X obtained in a second 
application of either of the two algorithms, then x ss = . . and correspondingly, if 



e^e^, where denotes the i-th column of I k , 



then clearly the diagonal entries of the matrices X and n T XIl will will be positioned 
in the same way along the diagonal. 
The fundamental observation is that 

(4.7) U T E st U = (^e^ej^j e s ej ^e^e^j = e 7r(s) e^ (t) = E n{sU{t) . 

This accounts for the permutations. Transposition is a little more complicated. The 
point is that after the diagonal variables in X are positioned, say Xi s is the ss entry of 
X for s = 1, . . . , k, as above, and if s ^ t and L(xi s ) fl L(x it ) = {x u , x v }, then either 
x u is st entry and x v is the ts entry, or the other way around. Thus, if 





T 
e n(l) 


k 


(4.6) n = 




=E 




T 


j'=i 
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(1) If x u is assigned to the st position of X, then the polynomial containing terms 
of the form 

must be placed in the st position in the array P. But this means that either x w is 
in the sr position or the rt position of X for some r other than s or t, since the 
ss, tt and st positions are already occupied. (If it is in the rt position, then x z 
will be in the sr position, so there is no loss of generality in assuming that x w 
is in the sr position of X.) Therefore, there is no loss of generality in assuming 
thatXu, e L(x ia ) (~]L(x ir ). Consequently, 

if x u is put in the st position of X, x w will be in the sr position. 

(2) If x u is placed in the ts position, then the polynomial displayed above must be 
placed in the ts position of the array P. Consequently x w must be in either the 
tr position or the rs position. But since it is in L(xi s ) D L(xi r ), the only viable 
option is that it is in the rs position of X: 

if x u is put in the ts position of X, x w will be in the rs position. 

Thus, transposition of the two variables in the first step after the diagonals are fixed, 
moves X to X T . 

Now suppose that there are two sets of single letter monomials and 
fee™ , • ■ ■ , bx r - k , where the monomials ax ik and bxj k must occur in the same polynomial 

in p. In this case it is possible for X to have diagonal variables x ^ , • • • , Xi k and for X 
to have diagonal variables Xj x , • • • , Xj k . Therefore, we must show that if this happens 
that X is pt equivalent to Y and that X is pt equivalent to Y. 

If Xi x , ■ ■ • , x ik axe the diagonal variables for X, then x^ , ■ ■ ■ , Xj k must be the diag- 
onal variables of Y and Xi a and Xj s must occur in the same diagonal position of X and 
Y. Similarly, if Xj x , ■ ■ • , Xj k are the diagonal variables of X and x ix , ■ ■ ■ , x ik are the 
diagonal variables of Y, then Xj 3 and Xi s must occur in the same diagonal position in 
X and Y. As in (4.7), let II be the permutation matrix such that the diagonal entries 
of n T Xil and Y are positioned in the same way along the diagonal. It follows that the 
diagonal entries of il T YTI and X are positioned in the same way as well. 

When Alg.l or Alg.2 determines the sets L(x is ) for X and L(xj s ) for X for 1 < 
s < k, it does so by considering terms of the form ax™~ l x u and bx n ~ l x v . One readily 
sees that for each term of the form that occurs in a polynomial in V, there 

is a corresponding term bx n ~ x x v that occurs in the same polynomial. Therefore the 
fundamental observation is that for each term .), there is a corresponding 

term x v E L(xj s ) and furthermore, if x u is placed in the st position in X, Algorithm 
PolyPos will place x v in the st position in Y. Similarly, if x v is placed in the st 
position in X, this correspondence ensures that Algorithm PolyPos will place x u in the 
st position in Y. Thus the nc variables X and Y are pt equivalent to X and Y. 

□ 
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Now that we have shown that there is a strong relationship between any two pairs 
of matrices constructed by our algorithms, we want to exploit this relationship to con- 
struct nc polynomials for matrix pairs determined by our algorithms. In particular, if 
X, Y and X, Y are two pairs of matrices determined by our algorithms for a given 
family V and p(X, Y) is an nc representation of V, we would like to conclude that 
there exists an nc polynomial p such that p(X, Y) is also an nc representation of V. 
We shall see in Lemma 4.6 below that this is true. 

Lemma 4.6 supplements Theorem 1.5 and is formulated in terms of a pair of auxil- 
iary nc polynomials that are expressed in terms of the notation introduced in (4.5): 



(4.8) p t (X,Y) = J2c a ,pmp',a'(Y,X) and p(X, Y) = p(Y, X). 

Then it is clear that 

Pt (X T ,Y T )=p(X,Y) T . 

Lemma 4.6. Suppose that p 1 , ■ ■ ■ ,p k 2 is a collection of polynomials V with an nc 
representation p(X, Y) satisfying the conditions in Theorem 1.4. If DiagParl or Diag- 
Par2, ParPosX and PosY generate a pair of matrices X and Y, then there exists a nc 
polynomial p such thatp(X, Y) is an nc representation ofV. 

Proof. If p satisfies the conditions in Theorem 1.4, then the algorithms DiagParl or 
DiagPar2, ParPosX, PolyPos and PosY can be applied to construct X and Y . Since 
p(X, Y) is an nc representation of V, by Theorem 4.5 and equations (4.3), (4.5), and 
(4.8) the following must hold: 



(4.9) 



(1) 


X 


= n T xn, 


Y 


= U T YU, 


(2) 


X 


= n T x T n, 


Y 


= U T Y T U, 


(3) 


X 


= n T rn, 


Y 


= n T xn, 


(4) 


X 


= n T r T n, 


Y 


= n T x T n 



p(X, Y) is an nc rep. 

p t (X,Y) is an nc rep. 

p(X, Y) is an nc rep. 

p t (X, Y) is an nc rep. 



□ 



Using the uniqueness results developed in this section, we can now prove Theo- 
rem 1.5. 

4.4.1. Proof of Theorem 1.5. 

Proof. Given thatp(X, Y) and p(X, Y) satisfy the conditions in Theorem 1.4, Algo- 
rithms ParPosl or ParPos2, ParPosX, PolyPos, and PosY can successfully determine 
the pairs X, Y and X, Y. Therefore, we may apply Theorem 4.5 to obtain the desired 
result. □ 

4.5. Determining p(X, Y) given X, Y and the positions of the polynomials in V. 

Once the matrices X and Y and the positions of the polynomials in V are determined 
by the previous algorithms, it remains only to find an nc representation for V . This 
rests on the following elementary observation, which is formulated in terms of the 
notation introduced in (2.3): 
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Lemma 4.7. The monomial m a ^(X, Y) is a kxk array of polynomials in the variables 
x%, . . . ,Xfc2 of degree \a\ + \/3\. 

Proof. This is immediate from the rules of matrix multiplication. □ 
4.5.1. Algorithm NcCoef. Suppose that 

(4.10) V = {pi, ■ ■ ■ ,Pk 2 } is a family of k 2 polynomials 

of degree d or less in the commutative variables x\, ■ • ■ , x 2 k 2 an d that these variables 
are positioned in X and Y and the polynomials in V are positioned in an array 

PA(i) ••• P\(k) 

(4.11) P= : ■.. : 

. P\(k 2 -k+i) ■ ■ ■ P\(k 2 ) . 

as in (1 .3). Then, since any nc polynomial p(X, Y) in the variables X and Y of degree 
d can be expressed in terms of monomials m a> p(X, Y) as 

(4.12) P(X,Y)= Ca,prn aj p(X,Y), 

H+\P\<d 

the NcCoef algorithm reduces to solving the system of linear equations 

(4.13) p(X,Y)=F, 

for the unknown coefficients c a ^. This system of equations has a solution if and only 
if there exists an nc polynomial p such that p(X, Y) is an nc representation of V. In 
view of Lemma 4.7, the addition of monomials m a> p(X, Y) of degree higher than d in 
(4.12) does not effect the solvability of (4.13). 

Remark 4.8. Algorithm NcCoef determines whether or not an nc representation exists 
for a given pair ofnc variables X and Y. Klep and Vinnikov [10] have been working on 
various elegant abstract characterizations of those sets of X, Y, and P which admit an 
nc representation. However, their work [10] does not address the issue of implementing 
tests for these characterizations. 

The following example illustrates how Algorithm NcCoef can be applied to deter- 
mine p(X, Y) for a given family V. 



Example 4.9. Let 



Pi 


= ±x\ + 


Pi 


= 4X1X2 


P3 


= 4xiX 3 


Pi 


= 4x 2 x 3 



?4 + X 3 Xq + x 2 x 7 + Qx 6 x 7 + 2x 4 Xg + 6Xg. 



and suppose that 

(4.14) X=( X ' X2 V Y=( X5 X «) and P = ( Pl P2 

\ x 3 x A J \ x 7 x 8 J \p 3 p 4 

Our goal is to find an nc representation p or to refute its existence. 
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Discussion Since the given family of polynomials is homogeneous of degree two, the 
first step of Algorithm NcCoef is to form the polynomial 

aX 2 + bXY + cYX + dY 2 = p(X, Y), 

and then set p(X, Y) = P. This yields four relations, one for each entry: 

ax\ + 0x2X3 + 6x1X5 + CX1X5 + dx\ + CX3XQ + dxQXY = p\ 

ax\Xi + 0x2X4 + CX2X5 + 6x1X6 + CX4X6 + 0x5X6 + 6x2X3 + dxQXg = P2 

0x1X3 + 0x3X4 + 6x3X5 + CX1X7 + 6x4X7 + 0x5X7 + CX3X8 + 0x7X3 = P3 

0x2X3 + 0X4 + 6x3X6 + CX2X7 + dx$X7 + 6x4X3 + +CX4X8 + dx\ = P4. 

Upon matching the coefficients of the left hand side of the first row in the preceding 
array with those of p 1 , we readily obtain the list of linear equations 

a = 4, 6 + c = 2, c=l and d = 6, 

the unknowns must satisfy; hence a = 4, 6 = 1, c = 1, d = 6. It is then easily 
checked that this choice of coefficients works for the remaining three rows of the array. 
Therefore, for this particular choice of X and Y and positioning of the polynomials 

p(X, Y) = AX 2 + XY + YX + QY 2 . 

□ 

4.5.2. Homogeneous sorting and Implementation of NcCoef. In order to implement 
the NcCoef algorithm, it is convenient to first sort each polynomial pj in V as a sum 

d 

i=l 

of homogeneous polynomials (pj)[i\ of degree i in which (pj)u] is the sum of terms in 
Pj that are homogeneous of degree i and is taken equal to zero if there are no such 
terms. 
Let 

(4.15) V i = {{ Pl ) [i] ,...,{p k 2)u} 

and, if at least one of the polynomials in Vi is nonzero, try to find a homogeneous nc 
polynomial representation Pi(X, Y) of degree i for V{. We shall refer to this procedure 
as homogeneous sorting. 

An obvious consequence of Lemma 4.7 is: 

Lemma 4.10. A family V ofk 2 polynomials of degree < d in 2k 2 commuting variables 
admits an nc representation p(X, Y) if and only if Vi admits an nc representation 
Pi(X, Y) with the same X and Y for each 1 < i < d. 

The primary advantage of homogeneous sorting is apparent when implementing Al- 
gorithms DiagParl, DiagPar2, ParPosX, PolyPos, PosY and NcCoef. One applies 
these algorithms to the nonzero family V% for each i separately. Typically, to save 
computational cost, choose i as small as possible in order to minimize computations. 
Then, once X, Y and the position of the polynomials are determined, it is easy to 
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fill in the coefficients of the terms in Pi(X, Y) by comparison with the array of terms 
{(pi)\i\, ■ ■ ■ , (Pk 2 )\i\} for the remaining choices of i, one degree at a time, just as in 
Example 4.9. It is important that the same X and Y are used for each choice of i. Now 
we give a cautionary example. 

Example 4.11. Ifpi, . . . , Pa are as in Example 4.9, then the set of polynomials 

(4.16) qx = pi + (x\ + x 2 x 3 )xi + (x x x 2 + x 2 x 5 )x 3 

(4.17) q 2 = p 2 + {x\ + x 2 x 3 )x 2 + (xix 2 + x 2 x 5 )x 5 

(4.18) q 3 = p 3 + {x 3 Xi + x 5 x 3 )xi + (x 3 x 2 + x\)x 3 

(4.19) q 4 = p 4 + {x 3 xi + x b x 3 )x 2 + (x 3 x 2 + xf)x 5 

does not admit an nc representation even though the terms of degree two admit an nc 
representation and the terms of degree three admit an nc representation. 

Discussion The terms of degree two are exactly the polynomials Pi, ■ ■ ■ ,Pi considered 
in in Example 4.9 and hence either lead to the representation considered there, or to an 
equivalent representation that corresponds to conjugation of the matrices X, Y and the 
polynomial array matrix P by a 2 x 2 permutation matrix II (to obtain n T XLI, n T FLI 
and n T Pn in place of X, Y and P), or transposition or to an interchange of X and 
Y. But in all these shufflings, the variables {xi, x 2 , x 3 , x 4 } will belong to one of the 
matrices and the remaining variables {£5, x 6 , x 7 , x$} will belong to the other. 

The polynomials q 1 , . . . , g 4 , will not admit an nc representation because the added 
terms come from 

r 1 3 

Xi x 2 

_x 3 x 5 _ 

which involves a mixing of the variables in the matrices X and Y that are obtained by 
analyzing p a , . . . , p±. □ 

4.5.3. Effectiveness ofNcCoef. The final task in our analysis of Algorithm NcCoef is 
to show that it will either be successful for all pairs X, Y produced our algorithms or it 
will fail for all such pairs. The next proposition (4.12) characterizes the nc polynomials 
p such that p(X, Y) is an nc representation of V when X and Y are determined by 
Algorithms DiagParl or DiagPar2, ParPosX PolyPos and PosY. Proposition 4.13 then 
insures that a given family V will either have an nc representation p(X, Y) for all X, Y 
produced by our algorithms or will have no representation for any such X and Y . 

Proposition 4.12. Suppose that a given family V ofk 2 polynomials in 2k 2 commuting 
variables admits an nc representation p(X, Y) and that the algorithms DiagParl or 
DiagPar2, ParPosX PolyPos and PosY produce the matrices X and Y . Then p belongs 
to the set W ofnc polynomials that is specified in Theorem 1.4. 

Proof. The coefficient conditions defining W are the exact conditions required for the 
algorithms DiagParl or DiagPar2, ParPosX, PolyPos and PosY to be successful. Fur- 
thermore, these algorithms will only be successful if these conditions are satisfied. 
Therefore, the assumption that our algorithms generate X and Y implies that if p ex- 
ists, then p G W. □ 
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Proposition 4.13. Suppose that V is a family of k 2 polynomials in 2k 2 commuting 
variables that satisfies (4.10). Let X, Y and X,Y be distinct pairs of matrices deter- 
mined by separate applications ofDiagParl or DiagParl, ParPosX, PolyPos andPosY. 
Then Algorithm NcCoefwill successfully determine an nc representation p(X, Y) ofV 
if and only if it successfully determines an nc representation p(X, Y)ofV. 

Proof. If the NcCoef algorithm produces an nc representation p(X, Y) of V, then 
Proposition 4.12 implies that p belongs to the set W specified in Theorem 1.4. Then 
Lemma 4.6 implies that there exists an nc polynomial p e W such that p(X, Y) is an 
nc representation of V. The reverse direction follows by a similar argument. □ 

4.6. The Size and Cost of NcCoef. We now determine the cost required to imple- 
ment NcCoef. For now, we neglect the cost to form the linear systems associated with 
NcCoef and only focus on the cost of solving these systems. 

For given arrangements a and A of the variables xi, ■ ■ ■ , x 2 fc 2 m X and Y and the 
polynomials pi, ■ ■ ■ , p k 2 in V we apply algorithm NcCoef to obtain a system of equa- 
tions in the undetermined variables c a ^ as in (4.12) and (4.13). By applying the method 
of homogeneous sorting as described in Section 4.5.2, we obtain (d — 1) systems of 
equations Mi in the unknowns c a p formed by equating 

(4.20) F\=pi(X,Y), 

where P A and p l a (X, Y) are formed from (4.11) and (4.12). 
To determine the cost of solving this linear system let 

(4.21) tij = number of monomials in P^L\, 

and set Tj = 5^i=i Uj. The number of noncommutative monomials in the variables 
X and Y of degree i is 2\ Therefore the number of unknowns c a p in our system 
corresponding to homogeneous terms of degree i will be T . It follows that 

Mi is a Ti x 2 l system of equations 

in the unknowns c Q/3 satisfying \a + @\ = i. 

When n > 2\ Mi will be overdetermined. The cost of solving this system using an 
LU decomposition is 

gi 

(4.22) fc 2 Ti4' arithmetic operations (see [7] §3.2.11) 

3 

We must solve such a system for each 2 < i < d, so the total cost to solve the linear 
system when each k 2 Ti > 2 % is 

d ( 8* \ 

(4.23) TotLinC CT A < ^ f Tj4 l — — J arithmetic operations. 

Similarly, when n < 2 l we get 

d 2 3i+1 

(4.24) TotLinC CT A < arithmetic operations (see [7] §3.2.9). 
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4.7. Final Results for One letter Algorithms. To this point we have developed the 
following algorithms: 

DiagParl: Partitions the commutative variables between the matrices X and Y and works 
under the assumption that <p(n, 0) ^ for some n > 2 and nip(n, 0) ^ 
ipin - 1,1). 

DiagPar2: Partitions the commutative variables between the matrices X and Y and works 
under the assumption that cp(n, 0) ^ for some n > 2 and (p(n, 0) ^ 
<p(n - 1, 1; Y) or <p(n, 0) ^ <p(Y; n-1, 1). 

ParPosX: Positions the variables in the matrix X if (p(n, 0) ^ for some n > 2 and the 
commutative variables are partitioned. 
PosPol: Positions the polynomials in the family if cp(n, 0) ^ for some n > 2 and the 
commutative variables are partitioned and positioned in X. 
PosY Positions the commutative variables in Y if (p(Q, n) ^ for some n > 1, the 
variables are partitioned and positioned in X and the polynomials are posi- 
tioned. 

NcCoef Given X, Y and the positioning of the pj in a matrix, this algorithm determines 
whether or not an nc p exists such that p(X, Y) generates the matrix containing 
the pj. 

Alg.l will refer to the sequence DiagParl, ParPosX, PosPol, PosY and NcCoef. 
Alg.2 will refer to the sequence DiagPar2, ParPosX, PosPol, PosY. and NcCoef. 



The next theorem summarizes the applicability of these two algorithms. 



Let iVC(4.25) denote the class of all nc polynomials p(X, Y) of degree d > 1 for 
which there exists integers s > , t > 1 and n > 2 such that 

f <p(s, t) ± 0, (p{n, 0) ^ 0, <p(n, 0) ^ ^(0, n), 
1 and np(n, 0) ^ p(n — 1,1). 

Theorem 4.14. Let V be a family ofk 2 polynomials in 2k 2 commuting variables. Then 
Alg.l yields an nc representation p(X, Y) ofV if and only if the given family V admits 
an nc representation in NC^s)- 

Proof. If Alg. 1 determines an nc representation p(X, Y), then Proposition 4.12 implies 
that p G iVC(4.25). Now suppose that V admits an nc representation in NCu 2 5)- Then 
Theorem 3.1, Theorem 4.1, Proposition 4.3, and Proposition 4.4 imply that Alg.l will 
successfully determine a pair of nc variables X and Y. Moreover, Lemma 4.6 implies 
that there exists an nc polynomial p such that p(X, Y) is an nc representation of V. But 
this implies that for this choice of X and Y that Algorithm NcCoef will be success- 
ful. Furthermore, Proposition 4.13 implies that for any pair of nc variables produced 
by Alg.l, algorithm NcCoef will be successful. Therefore, Alg.l will successfully 
determine an nc representation of V. □ 
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Let iVC(4.26) denote the class of all nc polynomials p(X, Y) of degree d > 1 for 
which there exists integers s > 0, t > 1 and n > 2 such that 



(4.26) 



<p(s,t) 7^ 0, <£>(n, 0) 7^ 0, y?(^, 0) 7^ 93(0, rz) , and either 
tp(n,0)^<p(n-l,l;Y), or p(n, 0) ^ <p(Y; n - 1, 1). 



Theorem 4.15. Let V be a family polynomials in 2k 2 commuting variables. ThenAlg.2 
yields an nc representation p(X, Y) ofV if and only if the given family V admits an 
nc representation in NC^ 2 q). 

Proof. If Alg.2 determines an nc representation p(X, Y), then Proposition 4.12 implies 
that p e iVC(4.26)- Now suppose that V admits an nc representation in NCu 26 y Then 
Theorem 3.2, Theorem 4.1, Proposition 4.3, and Proposition 4.4 imply that Alg.2 will 
successfully determine X and Y. Furthermore, Lemma 4.6 implies that there exists 
an nc polynomial p such that p(X, Y) is an nc representation of V. But this implies 
that for this choice of X and Y that Algorithm NcCoef will be successful. Further- 
more, Proposition 4.13 implies that for any pair of nc variables produced by Alg.2, 
algorithm NcCoef will be successful. Therefore, Alg.2 will successfully determine an 
nc representation of V. □ 

Remark 4.16. There are analogues of Theorems 4.14 and 4.15 in which NCu 25) is 
replaced by the class NCu 2 7) of all nc polynomials of degree d > 1 for which there 
exists integers s > 1, t > and n > 2 such that 



(4.27) 



<p(s,t)^0, f(0,n)^0, ^(n,0)7^(0,n) 
and n<f(0,n)^<f(l,n — l). 



and iVC(4.26) is replaced by the class NC^s) of all nc polynomials of degree d > 1 
for which there exists integers s > 1, t > and n > 2 such that 



(4.28) 



<f(s,t)^0, (p(0,n)^0, ip(n, 0) 7^ <^(0, n), and either 
(p{0,n)^(p{l,n-l]Y), or <p{0, n) ± <p(Y; 1, n - 1). 

5. Examples 



This section is devoted to a number of examples to illustrate the algorithms that were 
developed in earlier sections, as well as some variations thereof. 

Example 5.1. The set of A polynomials 

Pi = xix 5 + 5x1 + x 2 x 7 + 5x e x 7 
p 2 = XiXq + 52 5 a; 6 + x 2 x 8 + 5x 6 x 8 
p 3 = X3X5 + X4X7 + 5X5X7 + 5X7X 8 

Pa = X3X6 + 5x6X7 + X4X8 + 5xg 

in the commutative variables X\,...,xs can be identified with the entries in an nc 
polynomial p(X, Y) for appropriate choices of the 2x2 matrices X and Y. The 
objective is to find such an identification. 
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Discussion Since these 4 polynomials are homogeneous of degree 2, it is reasonable 
to look for an nc representation must be of the form 

p(X, Y) = aX 2 + bXY + cYX + dY 2 

for some choice of a, b,c,dE M. Moreover, since there are only two one letter mono- 
mials in the given family of polynomials: 

5x 2 inpi and 5xginp 4 , 

the corresponding variables must sit on the diagonal of either X or Y. We shall arbi- 
trarily place x 8 in the 11 position of X and x 5 in the 22 position of X. This in turn 
forces pi to be in the 22 position of p(X, Y) and p<± to be in the 1 1 position of p(X, Y) 
and forces a = 5 and d = 0. Furthermore, if A = X and B = Y with x 5 = x 8 = a, 
one of the other variables X{ = (3 and the remaining five variables equal to zero, then 
AB = BA and 

p(A,B) = 5A 2 + {b + c)AB, 

i.e., in terms of the notation introduced in subsection 1.3.2 with Cy = ip(i,j) for short, 
C20 — 5, Cn = b + c and c 02 = 0. Consequently, there are six possibilities: 
1) If Xi is the 12 entry of X, then 



a (3 


, B = 


"0 


0" 


a 









A 



2) If Xi is the 12 entry of Y, then 



and p(A, B) 



ha 2 lOa0 
5a 2 



A = 


'a 0" 

a 


, B = 


"0 



and p(A, B) = 


5a 2 Cnaf3 
5a 2 


3) If Xi is the 21 entry of X, then 








A = 


'a 0" 

(3 a 


, B = 


0" 



and p(A, B) = 


"5a 2 " 
10a/3 5a 2 



4) If Xi is the 21 entry of Y, then 



.4 




a 





, B 















a 




p 





5) If Xi is the 11 entry of Y, then 






A = 


a: 



0" 

a 




B = 


'ft 



0" 






6) If Xi is the 22 entry of Y, then 






A = 


a 


0" 




B = 


"0 


0" 







a 


? 





/3 





and p(A, B) 



and p(A, B) 



and p(A, B) 



'5a 2 

Cnaf3 5a 2 



5a 2 + cna/3 
5a 2 



5a 2 
5a 2 + cna/3 



Thus, if: 

(1) Xi is an off-diagonal entry of X, 2 polynomials will equal 5a 2 , one will equal 
10a (3 and one will equal 0. 

(2) Xi is an off-diagonal entry of Y, 2 polynomials will equal 5a 2 , one will equal 
cna(3 and one will equal 0. 
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(3) Xi is a diagonal entry of Y, one polynomial will equal 5a 2 , one will equal 
5a 2 + Cna/3 and two will equal 0. 

Correspondingly, if x 5 = x$ = a and Xi = (3 and all the other variables are set equal 
to zero, then the polynomials pi, . . . , p 4 assume the values shown in the following array 





i = 1 


% = 2 


i = 3 


i = 4 


i — 6 


% = 7 


Pi = 


5a 2 + af5 


5a 2 


5a 2 


5a 2 


5a 2 


ha 2 


P2 = 





a(3 








10a/3 





P3 = 








af3 








10ap 


Pi = 


5a 2 


5a 2 


ha 2 


5a 2 + af3 


5« 2 


ha 2 



Upon comparing the values of the polynomials p 1; . . . ,p 4 for the 6 possible choices 
of Xi = (3 with the possibilities (l)-(3) indicated just above, it is readily seen (from 
setting (3)) that X\ and X4 are diagonal entries in Y and c\\ — 1. (In fact since pi is in 
the 22 position of p(X, Y), x\ must be in the 22 position of Y, which forces X4 to be 
in the 11 position of Y .) The variables x^ and 23 will be off-diagonal entries of Y by 
setting (2); and hence the remaining variables xq and xj must be off-diagonal entries 
in X. A more detailed analysis would serve to position these last 4 variables in X and 
Y, after fixing one of them; see Remark 5.2. 

Remark 5.2. As p\ = X1X5 + ■ ■ ■ , p\ is in the 22 position of p(X, Y) and 25 is in 
the 22 position of X, it follows that x\ is in the 22 position ofY. Similarly, since 
Pi = X4X8 + ■ ■ ■ , Pa is in the 11 position ofp(X, Y) and x$ is in the 11 position of X, 
it follows that x 4 is in the 11 position ofY. Moreover, since the term 5x 6 x 7 in p\ (and 
P4) can only come from 5X 2 , it follows that Xq and x 7 must belong to X. Consequently, 
the remaining two variables, X2 and x$, must belong to Y. 

To go further, assume that x Q is in the 12 position of X. Then x 7 must be in the 21 
position of X and the polynomial 

P2 = x\Xq + ■ ■ ■ is in the 12 position ofp(X, Y). 

Therefore, X2 is in the 12 position ofY and the remaining variable £4 is in the 21 
position ofY. 

To better illustrate the algorithms, we shall return to the case where it is only known 
that xj G Y for j = 1, . . . , 4, Xj G X for j = 5, . . . , 8, x 8 is the 1 1 position of X and 
£5 is in the 22 position of X. Then L\ = {xq,x 7 } = L 2 and the positions of these 
two variables in X are not uniquely determined. We shall arbitrarily place x 6 in the 12 
position of X. Then x 7 must be in the 21 position and 

x\ + XqX 7 X 8 Xq + XqX 5 
X 7 Xg + £5X7 x\ + X 7 Xq 

Thus, P2 = 5xgXe + 5xqx$ + ■ ■ ■ must sit in the 12 position of p(X, Y). Therefore, 
x 2 x 8 is also in the 12 position of p(X, Y), which forces x 2 to be in the 12 position of 
Y. Therefore, 

„ _ Xg Xq 
X 7 £5 



5X 2 = 5 


x$ 


Xq 


= 5 




x 7 


X 5 





and Y 



X4 x 2 
x 3 x 1 
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It is now readily checked that 



Note! 



Example 5.3. The polynomials 

„2 



Pa P2 

_P3 Pi 

Pi P2 

P3 Pi 



5X 2 + XY. 



^ bX 2 + YX. 



x 1 + x 2 x 3 - x 3 x G + x 2 x 7 

XiX 2 + X 2 X A - X 2 X 5 + XiXq - 2 4 2 6 + X 2 Xg 

XiXj + X4X7 — X 3 X$ 



XiX 3 + X3X4 + X 3 X 5 
„2 



Pi = 

P2 = 
P3 = 

Pi = x 2 x 3 + x\ + x 3 x 6 - X 2 X 7 
admit an nc representation. 

Discussion Since these 4 polynomials are homogeneous of degree 2, it is reasonable 
to look for an nc representation must be of the form 

p{X, Y) = aX 2 + bXY + cYX + dY 2 

for some choice of a, b,c,d £ R, just as in Example 5.1. Moreover, since there are 
only two one letter monomials 



X'i 



inpx and x\ inp 4 , 



the corresponding variables are placed on the diagonal of X. Thus, a = 1 and d = 
0, and the substitutions X = A and Y = B in the candidate p(X, Y) for the nc 
representation yields the formula 

p(A, B) = c 20 A 2 + c n AB + c 02 B 2 , 

with c 2 o = 1, Cn = b + c and c 02 = 0. 

We shall arbitrarily place x 4 in the 11 position of X and x\ in the 22 position of X. 
This forces p 4 and p\ to be in the 11 and 22 positions of p(X, Y), respectively. Then, 
setting xi = x 4 = a and one of the other variables x { = (5 leads to the following sets 
of values: 

f 





x 2 = (3 


x 3 = (3 


x 5 = (3 


x& = (3 


x 7 = (3 


^8 = 


Pi 


a 2 


a 2 


a 2 


a 2 


a 2 


a 2 


P2 


2a(3 

















Ps 





2a(3 














Pi 


a 2 


a 2 


a 2 


a 2 


a 2 


a 2 



Since c 20 = 1 and c Q2 = 0, the discussion in Section 3.3 predicts 

2 polynomials equal to a 2 , one equal to 2a(3 and one equal to 0, or 
2 polynomials equal to a 2 , one equal to Cna/3 and one equal to 0, or 
1 polynomial equal to a 2 , one equal to a 2 + cna(3 and one equal to 0, 
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according to whether Xi is an off-diagonal entry of X, X{ is an off-diagonal entry of 
Y or xi is a diagonal entry of Y. Thus, c n = 0, x 2 and x 3 must belong to X, and 
x 5 , . . . ,x 8 belong to Y. 

Next, we shall arbitrarily place x 2 in the 12 spot of X. Then x 3 must be in the 21 
spot of X, and p 2 is in the 12 spot of p(X, Y). Now, having placed the entries in X 
and positioned the polynomials p 1: . . . , p 4 in p(X, F), it is readily seen that 



XiXq a term in p 2 > x i in the 22 position of X - 
X2X7 a term in p±, x 2 in the 12 position of X 



x 6 in the 12 position of Y 
> x-j in the 21 position of Y . 



However, it is not possible to fix the positions of x 5 and x 8 within Y from the available 
information. Thus, to this point we know that 



X 



X4 x 2 

X3 X\ 



and that either Y = 



X 5 Xq 
X 7 x$ 



or Y = 



X 8 X 6 
X-j X5 



Pi 


P2 


= x 2 + x 


X 8 Xq 




Xg x 6 


_P3 


Pl_ 


x 7 x 5 _ 




X-j x 5 _ 



and, since c 20 = 1 and cn = c 2 = 0, that there should be an nc polynomial of the 
form 

p(X, Y)=X 2 + a(XY - YX) 

for some a e K Since the coefficients of all the terms in p 1: . . . ,p 4 are ±1, it follows 
that a = ±1. It is then readily checked that 

X, 

X-j x 5 _ 

i.e., the second choice of Y works, the first does not. 

Example 5.4. Let k = 5 and n = 3 and suppose that xi is in the ii position for 
i = 1, . . . , 5 and that 



X 



The objective is to analyze the family of 25 commutative polynomials in 25 variables 
determined by X 3 to recreate X. As in the previous examples, the family of polyno- 
mials is the given data; we only include X here because we will not write out the full 
family determined by X 3 due to its prohibitive size. The reader is encouraged to create 
the family generated by X 3 using Mathematica and then follow along in the analysis 
to recreate X. 

Discussion If X is of the given form, then X 3 generates a 5 x 5 array of homogeneous 
polynomials of degree three, p 1: . . . ,p 2 5- Five of these polynomials will each contain 
exactly one term of the form x 3 . To simplify the exposition, we shall assume that the 
variables are indexed so that Xi is in the ii position of X for i = 1, . . . , 5. This in turn 
forces the polynomial that contains xf to be in the ii position of X 3 . The rest of the 
construction is broken into steps. 



x 1 


x 6 


x 8 


x w 


X12 


X-j 


x 2 


x 14 


xw 


xis 


Xq 


Zl5 


%3 


x 2 o 


X22 


X U 


x 17 




X4 


X24 


Xl3 


X\% 


X23 


X25 


x 5 
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1. Find the entries in Lj = Ri U Cj for % — 1, . . . , 5 by considering the terms xfxj 
that appear in X 3 , % — 1, . . . , 5 and j — 6, . . . , 25. For the given X we will obtain: 

L\ = {xq, Xj, Xs, Xg, X\o, In, £12, ^13}? 

L 2 = {:r 6 ,£ 7 ,£i4,£i5,£i 6 ,a;i 7 ,xi 8 ,xi9}, 

-^3 — {^8; ^9; ^14, ^15; ^20; #21 5 X 22 , ^23}) 

= {^10, ^11; %16, x 17, x 20, x 21, x 24, X25} 
L5 — { x 12, %13, x 18, x 19, x 22i ^23; ^24, x 2b] 

2. Position the entries in L\. 

Observe that L\ n L 2 = {x 6 ,a; 7 }. This means that one of these variables is in 
the 12 position and the other is in the 21 position. We shall assume that xq is in 
the 12 position and shall deduce the positions of all the other variables from the 25 
polynomials pi, . . . ,p 2 5, corresponding to X 3 . In particular, the assumption that x 6 is 
in the 12 position implies that the polynomial 

x\x& + x 6 xl + rri£ 6 £ 2 + XiX S Xi 5 + x i x w x n + x i x n x i9 H 

is in the 12 position of X 3 . But the remaining variables x\x s x t in that polynomial with 
off-diagonal variables x s and x t , s ^ t, will be in the 12 position of X 3 if and only if 
they are positioned in one of the ways indicated in the following array: 

x s 13 14 15 32 42 52 
x t 32 42 52 13 14 15 ' 

which is to be read as: 

either x s is in the 13 position and x t is in the 32 position, or vice versa 

The pair x s = x 8 and x t = x± 5 are subject to these constraints. On the other hand, 
since x 8 G L\ n L 3 it can only be in either the 13 position or the 31 position, whereas 
£15 G L 2 fl L 3 and hence can only be in the 23 position or the 32 position. Thus, the 
only viable solution to all these constraints is that 

x 8 is in the 13 position and x 15 is in the 32 position. 

Similarly, since xix 10 x 17 and x±xi 2 x w are in p 12 , whereas x w G LiflL 4 , x 17 G L 2 flL 4 , 
x i2 G L\ fl L 5 and x 19 G L 2 fl L 5 , it is readily seen that 

xio is in the 14 position and x 17 is in the 42 position 

and 

X12 is in the 15 position and x w is in the 52 position. 

Moreover, now that x 6: x%, x w and x\ 2 are positioned, we consider the following mono- 
mials in p 12 : 

xqXsXq =>■ Xg is in the 31 position of X 
xqXiqXu =>- Xu is in the 41 position of X 
xqXi 2 xi 3 =>- x 13 is in the 51 position of X. 
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Remark 5.5. The preceding calculations exploit the entries in the polynomial in the 
12 position ofp(X, Y) to calculate terms in R lt C\ and R 2 . A variant of this is to just 
fill in Ri and then to rely on successive steps to fill in R 2 , . . . , R5, one row at a time. 
At the other extreme, it is also possible to use all the entries in this polynomial to fill in 
the whole matrix; see item 6, below. 

3. Position the remaining entries in L 2 by the algorithm. 
To this point we know that 



X\ 


Xq Xg 


#10 


%12 


x 7 


x 2 






x 9 


%15 x 3 






Xll 


Xyj ■ 


X4 




Xl3 


X 19 ■ 




x 5 



and it remains to fill in the dots. Since the polynomial 

x 7 x\ + xlx 7 + X 2 XuX 9 + X 2 XiqX U + X 2 Xi S X 13 H 

sits in the 21 position of X 3 and the positions of x 9 , x n and x 13 are known, it follows 
that 

X14 is in the 23 position, x±e is in the 24 and xi 8 in the 25. 
Moreover, since 

{x u , ^15} = L 2 fl L 3 and x 14 is in the 23 position 

it follows that x\$ is in the 32 position. Similarly, since 

{^16, xn} = L 2 f] L A and x w is in the 24 position 

it follows that x 17 is in the 42 position. Much the same argument based on the position 
of #i8 and the observation that L 2 flL 5 = {x%s, x± 9 } shows that x±g is in the 52 position 
of X. This completes the positioning of the variables in L 2 . 

4. Position the remaining entries in R 3 using Algorithm ParPosX outlined in Sec- 
tion 4.1 

Since the polynomial 

£3X9 + X 9 xj + £3X15X7 + X 3 X 20 Xn + X 3 X 22 X!3 H 

is in the 31 position of X 3 and the positions of x 3 , X7, xn and x i3 are known, it is 
readily checked that x 20 is in the 34 position and x 22 is in the 35 position. 

5. Position the remaining entries in X using Algorithm ParPosX. 
Since the polynomial 

x\x\\ + X\\x\ -\- X4X17X7 + X4X21X9 + X4X24X13 -(-••• 

is in the 41 position of X 3 and the positions of and X13 are known, it follows 

that x 2i is in the 43 position and x 24 is in the 45 position. Similarly, as the polynomial 

X5X13 + X13X? + X5X19X7 + x 5 x 23 a;9 + ^5^25^11 H 
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is in the 51 position of X 3 , the positions of x$, xj, xg and x±i serve to locate x 19 in 
the 52 position, x 2 3 in the 53 position and x 25 in the 54 position. This completes the 
computation via Algorithm ParPosX. 



6. Variations on the theme. 

The preceding steps fill in the entries of the partially specified matrix X in (5.1) by 
studying the entries in 5 of the given 25 polynomials by sweeping along rows. It is 
also possible to fill in X by using all the entries in the polynomial 



q 



xlx 6 + x 6 xj + xix e x 2 + xix 8 xi 5 + xix w x 17 + xix 12 x 19 + 



that sits in the 12 position of X 3 : 

x 6 x 8 Xg a term in q, x 6 in the 12 spot, x 8 in the 13 spot = 

x 6 x w xn a term in q, x e in the 12 spot, x 10 in the 14 spot = 

XGX12X13 a term in q, x 6 in the 12 spot, X12 in the 15 spot = 

XqXuXi 5 a term in q, x e in the 12 spot, x 15 in the 32 spot = 

xqXiqXi7 a term in q, x 6 in the 12 spot, xn in the 42 spot = 
a term in q, x 6 in the 12 spot, xig in the 52 spot = 

x 8 x 2 o^i7 a term in q, x 8 in the 13 spot, x 17 in the 42 spot = 

£8X22X19 a term in q, x$ in the 13 spot, x\g in the 52 spot = 

£10X15X21 a term m <?> x w hi the 14 spot, X15 in the 32 spot 

X12X15X23 a term in q, x\ 2 in the 15 spot, x i5 in the 32 spot 

£12X17225 a term in Q, X12 in the 15 spot, xn in the 42 spot 



x 9 in the 31 spot 

> £n in the 41 spot 

> X13 in the 51 spot 

> Xu in the 23 spot 

> xi 6 in the 24 spot 

> Xis in the 25 spot 

> x 2 o in the 34 spot 

> X22 in the 35 spot 
=>■ x 2 i in the 43 spot 
=>■ x 2 3 in the 53 spot 

x 2 5 in the 54 spot 



Therefore, the single remaining variable x 2 4 must be in the 45 spot. 

Remark 5.6. If, in Step 2, x 6 is placed in the 21 position and Xj is placed in the 12 
position, then the algorithm will generate the transpose X T of the matrix X that is 
specified in the statement of Example 5.4. 



6. Cost of Alg.2 vs. a Brute Force Approach 

Now that we have developed our one-letter methods for determining an nc represen- 
tation and determined which families they will work for, we would like to get some 
idea of their cost. As a sample we get some rough cost estimates for Alg.2 and omit es- 
timates for Alg.l. As Alg.l and Alg.2 differ only in the partitioning algorithms used, 
we suspect that the costs for Alg.l will be similar to the costs for Alg.2. The goal 
of determining these estimates is to compare the gains in efficiency that our methods 
provide over the more direct, brute force approach described in §1.6. As usual, we 
assume V is a family containing k 2 polynomials of degree d or less in the commutative 
variables xi, ■ ■ ■ , x 2 k 2 - 
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6.1. Brute Force. Recall that in Section 4.6 we determined the cost of solving the 
linear systems determined by NcCoef. As the Brute Force Method merely applies al- 
gorithm NcCoef for each arrangement of the commutative variables and polynomials, 
our work is nearly finished. However, recall that in our cost analysis of algorithm Nc- 
Coef we omitted the cost to form the systems. For the sake of completeness, we briefly 
mention how to form these linear systems using a function like the one in Mathematica 
called CoefficientList[] to get an idea of the cost. 

The function CoefficientList[] in Mathematica takes as an input a polynomial and 
a collection of commutative variables and outputs the coefficients of the polynomial 
associated with the given variables. For example, the command 

A[i\ = CoefficientList[£>j, {xi, • • • , x 2k i}] 

creates a multidimensional array A[i\, where the entry A[z](ii, • • • , i 2 k 2 ) corresponds 
to the coefficient of the monomial x^ ■ ■ ■ x^jS in Pi. The starting point for our imple- 
mentations of both the brute force approach and Alg.2 will be to call CoefficientList[ 
for each polynomial in the given family V. We note that any implementation of our 
algorithms will require a function similar to CoefficientList[] in order to of access the 
coefficients of the family pi, • • • ,p 2 k 2 in a systematic way. 
Let Mi be the system formed by equating 

(6.1) F\=pi{X,Y), 



where P A and p % a [X, Y) are formed from (4.11) and (4.12) respectively, by homoge- 
neous sorting. Then letpAQ) be an element of the matrix P A , and let qj be an element of 
Prr(X, Y) with undetermined coefficients c a p. Letting A[j] = CoefficientList[gj, {x±, ■ ■ ■ , £2fc 2 }] 
for each 1 < j < k 2 , and B[j] = CoefficientList[pA(j), {xi, ■ ■ ■ , x 2 k 2 }], we obtain the 
system associated with Eq. (6.1) by setting 

(6.2) A\j](i h ■■■ , i 2k i) = B[j](z u • • • , t 2k 2) (1 < j < k 2 ) 

for (i h ■ ■ ■ , i 2k 2) satisfying J2l=i ^ = *• 

Excluding the cost of iterating through the terms in Eq. (6.2), the cost to form the 
systems of equations associated with NcCoef involves k 2 calls to the function Coeffi- 
cientList[]. In addition to this, recall that the cost to solve the linear systems formed by 
NcCoef for a given arrangement a of the variables and A of the polynomials satisfies 

d ( 8 i \ 

(6.3) CostNcCoef^A < ( ^4* J arithmetic operations, 

provided > 2*. The formula is similar if r < 2*. 

Recall that there are (2k 2 ) \ arrangements a of x±, ■ ■ ■ , x 2 fc 2 in X and Y and (k 2 ) \ ar- 
rangements A of pi, ■ ■ ■ , Pk2. We note that we must call CoefficientList[] once for each 
Pj in the given family V and once for each undetermined polynomial qj determined by 
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a. Therefore if r > 2\ the cost satisfies 

(6.4) CostBruteForce < (2k 2 )\(k 2 )\ ~~ \jj ^hmetic operations 

+ ((A; 2 + (2A; 2 )!) Calls to CoefficientList[]) . 

6.2. Cost of Algorithm 2. In this section we roughly analyze the cost associated with 
each subroutine required to execute Alg.2. Then in Section 6.2.7 we summarize these 
costs and compare them to the brute force cost determined in Eq (6.4). 

6.2.1. Algorithm DiagA: Determining the Diagonal Entries. The first step in Alg.2 is 
to determine the commutative variables that occur in one-letter monomials of the form 
axf . Here we give a rough outline of a procedure based on Lemma 2.7 that will do 
this, and then analyze the cost. 

As in the case of the brute force algorithm and algorithm NcCoef, we implement 
this algorithm using the function CoefficientList[]. Again, we let 

A[i\ = CoefficientListfpj, {x 1 , ■ ■ ■ , x 2 k 2 }]- 

Using the coefficient matrix A[i] for each 1 < i < k 2 , we inspect the coefficients of the 
single letter monomials of degree n, where 2 < n < d. For example, the coefficients 
of the one letter monomials ax r [ and bx% in pi are given by 

A[i](n,0,--- ,0) = a and A[i](0, n, 0, • • • , 0) = b. 



The core of Algorithm DiagA is to: 

Fix a degree 2 < n < d. Then iterate through all of the coefficients of one letter 
monomials of degree n in V. That is, iterate through the expressions 

(6.5) A[l](n,0... ,0), ••• ,A[1](0,--- ,0,n) 

A[2](n,0,... ,0),... ,A[2}(0,--- ,0,n) 

A[k 2 ](n,0 r -- ,0),-.. ,A[k 2 }(0 r -- ,0,n) 

and verify that exactly 0, k or 2k of the above expressions are nonzero and that they 
assume at most two distinct, nonzero values. If this is satisfied, then store the nonzero 
coefficients and store the associated diagonal variable in one of two lists based on its 
coefficient. If this is not satisfied then quit. 

6.2.2. Cost Estimate for DiagA. The cost of algorithm DiagA is easy to bound and 
can be expressed in terms of the parameters k and d and the number of calls to Coeffi- 
cientList[]. Define an equality check to be a test if two given numbers are equal. 

Algorithm DiagA begins with k 2 calls to the function CoefficientList[], one for each 
polynomial p^. Iterating through the (2k 2 )(k 2 ) coefficients listed in (6.5) for a given 
degree 2 < n < d and performing a small number Ca of equality checks for each 
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coefficient requires no more than (d— l)(k 2 )(2CAk 2 ) = 2(7^(^ — 1) k A equality checks. 
Therefore, the total cost for applying DiagA is 

(6.6) TC = (C A (d - l)k 4 equality checks) + (k 2 Calls to CoefficientList[]) . 
Here Ca is certainly no greater than 5. 

6.2.3. Cost ofDiagParl. Here we determine an upper bound for the cost of Algorithm 
DiagPar2. We assume that DiagA has already been executed and that CoefficientList[] 
has been called for each polynomial pi E V. See §6.2.1 for a definition of these 
procedures. 

Algorithm DiagPar2 utilizes the coefficients of monomials of the form 
where x, ; is a diagonal variable and X{ ^ Xj, to partition the variables between X and 
Y. The bulk of the cost of this algorithm will lie in iterating through the coefficients of 
the two-letter monomials of this form. We fix a degree 2 < n < d, a polynomial pj and 
a diagonal element Xi (which can be determined using DiagA) and iterate through the 
coefficients of monomials inpj of the form x r l' 1 xi, where x\ is a non-diagonal element; 
therefore, we iterate through 2k 2 — 2 coefficients for each diagonal element. For each 
coefficient we perform no more than Cp ar 2, a small number, of equality checks to 
ensure that the conditions in Lemma 2.12 are satisfied and to determine whether xi is 
in X or Y, 

If at any point we find that these conditions are violated or we cannot find an n such 
that the two-letter monomials ax^~ l Xj satisfy the conditions consistent with (3.10), 
we terminate the iteration and conclude either that V has no nc representation or that 
the algorithm is inconclusive. As we are only considering k diagonal variables, k 2 
polynomials and d — 1 different degrees that need to be considered, we have that 

(6.7) TotCostDiagPar2 < Cp ar2 (d — 1) k 5 equality checks. 
Here C Par2 < 5. 



6.2.4. Cost of ParPosX. We now determine a rough bound on the cost of Algorithm 
ParPosX. Again, we assume Algorithm DiagA and DiagPar2 have been executed and 
that CoefficientList[] has been called for each of the polynomials pi e V. 

Algorithm ParPosX positions the variables in X by using lists L» defined in (4.1) 
and three letter monomials of the form x™~ 2 x;x m , where Xi is a diagonal and x\ ^ x m 
are non-diagonal elements in X. We observe that we can build the lists Lj when we 
iterate through the two letter monomials ax" _1 Xj in DiagPar2, so other than the cost of 
DiagA, DiagParl or DiagPar2, the primary cost associated with this algorithm lies in 
iterating through the three-letter monomials. For a given non-diagonal polynomial pj 
and diagonal variable Xi, we must iterate through the coefficients of [k 2 — k)(k 2 — k — 1) 
such monomials, and for each nonzero coefficient (of which there will be at most k — 1), 
iterate through a list of length 2k — 2 and perform a fixed number of equality checks 
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which we will designate by Cp ar p- Therefore, 
(6.8) 

TotCostPar2Pos =C ParP {k 2 - k)k[({k 2 - k) (k 2 - k - 1)) - k] + 

Cp ar p(k 2 — k)k{k — I) (2k — 2) equality checks 

< 3Cp ar pk 7 equality checks. 

Here Cp ar p < 10. 

6.2.5. Cost of PosPol. As in the case of Algorithm DiagPar2, the cost of Algorithm 
PosPol is effectively reduced to the cost of iterating through two-letter coefficients. 
As in DiagPar2, we inspect two-letter monomials of the form ax™ Xj, where Xi is a 
diagonal element of X and Xj is a non-diagonal element of X. For a fixed degree 2 < 
n < d, polynomial p m and diagonal element Xi of X, we iterate through the coefficients 
of all possible two-letter monomials in p m , which allows for k 2 — k possibilities. Once 
we find a coefficient a / corresponding to the monomial ax™~ x Xj, we immediately 
determine the position of p m using the position of xi and Xj. Therefore 

(6.9) TotalCostPolyPos < C PosP (d - l)k 5 equality checks, (C PosP < 10). 

We note that this is probably a gross overestimate of the cost of PosPol and that one 
does not need to iterate through all of the two-letter monomials again. It seems pos- 
sible that we could store additional data associated with the non-diagonal entries Xj 
when we iterate through the two-letter algorithms in DiagPar2 in order to position the 
polynomials. 

6.2.6. Cost ofPosY. We now determine a bound for the cost of PosY for families V E 
NC3.11. We again assume that all necessary arguments to implement PosY have been 
executed so that the commutative variables are positioned in X and the polynomials 
are positioned. Therefore we assume that half of the commutative variables have been 
positioned in X and that the polynomials pi, • • • ,Pfc2 have been assigned positions in 
a k x k matrix. The goal is to to position the remaining variables in Y 

The first step is to determine the diagonal variables of Y. PosY uses the fact that 
(p(s, t) 7^ for some s, t such that 2<s+t=n<d and t > 0. Then we look at two 
letter monomials of the form ip(s, t)xfxj in the diagonal polynomials pj, where Xi is a 
diagonal variable of X. This will allow us to conclude that xi is a diagonal variable in 
Y that must be located in the same position as Xi is in X. To implement this step, for 
a fixed degree n, diagonal polynomial pj and corresponding diagonal variable Xi in X, 
and s, t satisfying s + t = n,we must iterate through k 2 — k coefficients corresponding 
to two letter monomials of the form xfx\, where xi is in Y. Therefore 

(6.10) CostDiagY < C PosY (d - 1)0) (k 2 - k) equality checks 

< Cp OS Y(d 3 k 3 ) equality checks, 
where Cp os y dominates the operation per step count. 
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Finally, to determine the position of the off-diagonal variables of Y, we iterate 
through three-letter monomials of the form where X; and x t are diago- 

nal variables of X and Y respectively, and x v is a non-diagonal term of Y. We observe 
that if such a term occurs in a polynomial pj that is in the qr -position of the polynomial 
matrix and Xj and x\ are in the qq position of X and Y respectively, then x v must be in 
the qr-position of Y. Again, the bulk of the cost lies in iterating through these terms. 
For a fixed non-diagonal polynomial pj and fixed diagonal variables Xi and xi in the 
gg-position in X and Y, we consider the coefficients of k 2 — k such terms and fewer 
than Cp OS Y equality checks each . Therefore, 

(6.1 1) CostNonDiagY < C PosY {k 2 - k){k){k 2 - k) equality checks 

< Cp OS Y(k 5 ) equality checks. 

Thus, 

(6.12) TotalCostPosY < C PosY {(d 3 k 3 ) + (k 5 )) equality checks 
Here Cp os y < 10. 

6.2.7. Cost of Algorithm 2. The following table lists the subroutines that make up 
Alg.2, their cost in terms of the number of operations and calls to CoefficientList[], 
and the section in which the cost of the subroutines was determined. Recall that the 
parameter represents the total number of commutative monomials of degree 2 < i < 
d in all polynomials in V, and C DiagA-, Cp ar2 , Cp ar p, Cp os p, Cp os y are all constants 
bounded by 10. 



(6.13) 



Algorithm 


Operations 


Calls to CoefficientList[] 


DiagA 


C DiagA dk A 


k 2 


DiagPar2 


Cp ar 2 dk b 





ParPosX 


CparP k 1 





PosPol 


CposP dk b 





PosY 


C PosY (d*k* + P) 






Combining this with the cost of applying algorithm NcCoef described in §4.6, we 

obtain 

(6.14) 

CostAlg2 < 10 (k 7 + (2d + l)k 5 + dk 4 + d 3 k 3 ) equality checks 

d ( 8 i \ 

+ I Tj4* J arithmetic operations + k 2 Calls to CoefficientList[], 

i=2 ^ ' 

provided T- t > 2 % for each i. When Tj < 2 % for each i, we get 

(6.15) CostAlg2 < 10 (k 7 + (2d + l)k 5 + dk 4 + d 3 k 3 ) equality checks 

d 2 3i + 1 

+ arithmetic operations + k 2 Calls to CoefficientList[]. 
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6.2.8. Comparison of Costs. A comparison of the bounds (6. 4), (6. 15) and (6.14) shows 
the benefit of our Algorithms. We first observe that we need to use the function Co- 
efficientList[] (2k 2 )\ fewer times using Alg.2, which is a huge savings given that the 
cost to use CoefficientList[] could potentially be very expensive. Even if we neglect 
the cost of this function, we see that for large d and k, the cost to form and solve linear 
systems using NcCoef dominates both the cost for the Brute Force Method and Alg.2. 
By exploiting the structure of the polynomials in the given family V and performing 
on the order of 

k 7 + 3dk 5 + d 3 k 3 operations, 

we have effectively reduced the cost from solving {2k 2 )\{k 2 )\ such systems to a single 
system. This is a vast improvement. Also, to rule out the existence of an nc represen- 
tation using the Brute Force method we must check all of these cases and verify that 
they fail. Much to the contrary, Algorithm 2 is likely to determine non existence very 
early when applying it. 

7. Families that may not contain a term of the form ax" with n > 1 

Sections 3, 4 and 6 present algorithms for solving our nc representation problem 
when at least one of the given polynomials in V contains a term of the form ax™ with 
a G ]R\{0}andn > 1. This section treats V which contain no one letter monomials but 
which do contain terms of the form axfx*, a ^ 0. Recall Lemma 2.8 and Lemma 2.10 
dealt with patterns two letter monomials in an nc representable V must obey. 

The first step in developing these two letter algorithms is to determine the diagonal 
elements given the existence of two letter monomials. Lemmas 2.16 and 2.17 present 
conditions under which the presence of terms that are > -equivalent to xfx*- allows us 
to determine dyslexic diagonal pairs. The next step is to partition the k dyslexic pairs 
{x ii; Xj t }, ■ ■ ■ , {x ik , Xj k }, i.e., to determine which elements are on the diagonal of X 
and which are on the diagonal of Y. An application of Lemma 2. 19 is usually sufficient 
to partition these diagonal pairs. 

Once the diagonal variables are determined and partitioned, we can reduce the prob- 
lem to one that is manageable for the single variable algorithms by taking derivatives 
of the polynomials in V with respect to the diagonal variables. This process will be 
known as the SVR algorithm. The next section discusses this reduction. 

We shall not present a cost analysis of our two-letter methods in this section, since 
one can see as they unfold that they are clearly far superior to Brute Force. For one 
thing the core of our two-letter procedures are reductions to our single letter algorithms 
(such as Alg. 2 whose cost is vastly less than of Brute Force). 

7.1. Reduction to one letter algorithms: SVR Algorithm. In this section we de- 
velop the SVR (Single Variable Reduction) Algorithm, which can be used to reduce 
a family of polynomials containing terms that are > -equivalent to xfx* to a family of 
polynomials to which the single variable algorithms developed earlier apply. The next 
lemma contains the key idea that underlies this algorithm. 
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Lemma 7.1. Ifpi, . . . ,Pk? is a family of polynomials in the 2k 2 commuting variables 
Xi, . . . , X2fc2 that admits an nc representation p(X, Y) of degree d > 2 such that Xi r is 
in the ab position of X and x i± , . . . , x ik sit on the diagonal of X, then: 

(1) The polynomials 

dpi dp k 2 
dx ir ' " ' ' dx ir 

admit the nc representation 

(7.1) Jim j(p{X + tE ab , Y) - p(X, Y)). 

(2) The family of polynomials q m (xi, . . . , x^), m = 1, . . . , k 2 , defined by the for- 
mula 

q m (xi, • • • , x 2fc 2) = ^ J—Pm(xu ■ ■ ■ , X2k*), rn = 1, . . . , fc 2 , 

s=l ls 

admits the nc representation 

(7.2) lim \{p{X + tl k , Y) - p(X, Y)). 
Proof. The first assertion follows from the observation that 

]mi^{p{X + tE ab ,Y)-p{X,Y)) 
is equivalent to computing 

ry~\ I /y . nn , nn . I T rjn . nn . \ ff\ I ff . ff . ff . If . If 1 . \ 

y FmV x «i ) • • • ) ^V-l ' ,ht r ' 1 Jj lr+1 J ' ' ' 1 l 2k 2 ' t'm\ h %\ ; I • • • J ^Ir-l 5 1r ' ^V+l 1 ' • • J l 2k 2 ' 

t^o t 

Pm \p^i\i ■ ■ ■ i I ■ 



OXi, 



The second assertion follows from the first by a straightforward calculation. □ 

To ease future applications of Lemma 7.1 when x^, - ■ ■ , x ik are the diagonal ele- 
ments of X, it is convenient to introduce the notation 

k d 

(7.3) {T x Pm){x\, ■ ■ ■ ,x 2 k 2 ) = ^2 Q — Pm^i)'-' ,% 2 ), m = l,...,k 2 . 

s=l ls 

Remark 7.2. Repeated application of the formulas in Lemma 7.1 serves to reduce nc 
expressions in X and Y to expressions in the single variable Y. Thus for example, if 
Pi, . . . ,p k 2 is a family of polynomials corresponding to 

p(X,Y) = XY 2 X 2 , 

then the family of polynomials T 3 p m , m = 1, . . . , k 2 , corresponds to the polynomial 

3\p(I,Y) = 3\Y 2 . 
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In this way it is possible to eliminate the dependence on the k 2 variables in X by dif- 
ferentiating the given family of polynomials just with respect to the k diagonal entries 
ofX. 

Let D x ,i k p(X, Y) (resp., Dy t i k p{X, Y)) denote the directional derivative of p(X, Y) 
with respect to X(resp., Y) in the direction of the identity 1^. More generally, let 
D XIk p(X, Y) (resp., D YI p(X, Y)) denote the i-th directional derivative of p(X, Y) 
with respect to X (resp., Y) in the direction of Ik- Then Lemma 7.1 states that if the 
variables x^, ■ • • , Xi k are the diagonal elements of X (resp., Y), then Dx,i k p(X, Y) 
(resp., D Y j k p(X, Y)) is a representation of the family 

g n = T x pn (resp. g n = T Y p n ), 
where Tx (resp., Ty) is the operator defined in (7.3). 

Remark 7.3. By repeated application of Lemma 7.1, we have that D\ Ik p(X,Y) 
(resp., D l Y j k p(X, Y) ) is an nc representation of the family defined by 

g n = T l x p n (resp., g n = T Y p n ) 

Theorem 7.4 (SVR Algorithm). Let pi, ■ ■ ■ ,p k 2 be a homogeneous family of polyno- 
mials with an nc representation p(X, Y) of degree i + j where i > 2 and j > 2. 
Suppose that the k partitioned diagonal pairs 

are known and tp(i,j) ^ 0. Also assume that 

(7.4) jtp{i,j)^(i + l)tp{i + l,j-l) 
or 

(7.5) iip(i,j)^{j + l)(p{i-l,j + l). 

Then if (7.4) (resp., (7.5)) holds the variables can be partitioned by Algorithm Diag- 
Parl and then positioned in Y ( resp., X ) by Algorithm ParPosX applied to D x 7 p(X, Y) 
(resp.,D YJk p(X,Y)). 

Proof. Since the diagonal pairs are partitioned, we may assume that x ir , . . . , x ik are 
the diagonal entries of X and xj x , . . . , Xj k are the diagonal entries of Y. Moreover, the 
nc polynomial 

i+j 

p(X,Y) = J2ls(X,Y), 

s=0 

where q s (X, Y) denotes the sum of the terms in p(X, Y) that are of degree s in X and 
degree i + j — s in Y. Consequently, 

i+j 

(D x>Ik p)(X,Y) = ^(D x , Ik qs)(X,Y) 

(7.6) s= ° 

= <p(i,j)UY' + (D l XJk q i+1 )(X,Y) + lPxj k *e){X>Y). 

s=i+2 
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Only the first two terms on the right in the second line of (7.6) will contribute two 
letter monomials of the form ax 3 A~ x u . If X is replaced by A, Y is replaced by B and 
AB = BA, then these two terms will be equal to 

i\ if(i,j)B j + (i + 1)! <p(i + 1, j - l)AB i_1 

Therefore, if j(H(p(i,j)) ^ {i + 1) + 1, j — 1)> then we may apply DiagParl to 
partition the variables. But this is the same as (7.4). Moreover, in view of Theorem 

4.1, Algorithm ParPosX, will then serve to position the variables in Y. 

Similarly, if (7.5) is in force, then DiagParl and ParPosX applied to D d YI will serve 
to partition the variables and to position them in X. □ 

Remark 7.5. Theorem 7.4 remains valid if (7.4) and (7.5) are replaced by the condi- 
tions 

tp(i,j)^(i + l)tp(i + l,j-l;X) or <p(i,j)^(i + l)<p(X;i + l,j-l). 

and 

<p(i,j)^{j + l)<p(i-l,j + l:Y) or <p(i,j)J:(j + l)<p(y',i-l,j + l), 

respectively, but with DiagPar2 in place of DiagParl. (In the first case, DiagParl is 
applied to D l X I j>, in the second, it is applied to D YI p. 

The algorithm outlined in Theorem 7.4 will be called the SVR (Single Variable 
Reduction) Algorithm. 

7.2. Algorithms based on two letter words. In this section we consider algorithms 
for families of polynomials V that contain two letter words that are > -equivalent to 
xfx*j with % 7^ j and s > t > 2. It is convenient to separately analyze the three 
mutually exclusive cases 

s >t + 2, s = t + 1, and s = t. 

Our approach is to use either Alg.l or Alg.2 in combination with the SVR algorithm. 
Recall that Alg.l refers to the sequential application of DiagParl, ParPosX, PosPol and 
PosY, whereas Alg.2 refers to the sequential application of DiagPar2, ParPosX, PosPol 
and PosY. 

7.2.1. Two letter monomials \> -equivalent to x s u x l v with s > t + 2. This subsection 
contains two results that provide nc representations for families containing two letter 
monomials with s > t + 2 > 4; the first is based on Alg.l, the second on Alg.2. 

Let NC(j.7) denote the class of all nc polynomials p(X, Y) of degree d such that for 
some set of integers s, t with 

r s>t + 2>4, (f(s,t)^0, f(s,t) j£ tp(t,s) and either 

iLI) \ Ms, t)^(s + l)<p(s + 1, t - 1) or s<p(s, t)^{t + l)<p(s - 1, t + 1). 

Proposition 7.6. Let p±, . . . ,p^ be a family V of polynomials in Ik 2 commuting vari- 
ables. Then the SVR algorithm coupled with Alg.l will yield an nc representation of 
p(X, Y) ofV with p in NCft.Y) if and only if the given family V admits an nc repre- 
sentation in this set. 
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Proof. In view of the assumptions in the first line of (7.7), Lemma 2.16 may be ap- 
plied to obtain the partitioned diagonal pairs. The constraints in the second line of 
(7.7) then insure that SVR algorithm coupled with Alg.l serves to partition and po- 
sition the variables and to position the given set of polynomials. More precisely, if 
tip(s, t) (s+ l)(p(s + 1, t — 1) the differentiation is with respect to X; if sip(s, t) ^ 
(t + l)v?(s — 1, t + 1), then the differentiation is with respect to Y. Moreover, if the 
differentiation in the SVR algorithm is with respect to X, then the ParPosX Algorithm 
will position the variables in Y and the PosPol Algorithm will position the polynomi- 
als. On the other hand, if the differentiation in the SVR algorithm is with respect to Y , 
then the ParPosX Algorithm will position the variables in X. Finally, if the given fam- 
ily does not admit an nc representation, then the indicated algorithms cannot produce 
it. " □ 

Let NC(j,s) denote the class of all nc polynomials p(X, Y) of degree d such that for 
some set of integers s, t with 

s>t + 2>4, ip(s,t)^0, (p(s,t) ^ (p(t,s) and 
<p(s, t)^(s + l)<p(s + l,t-l;X) or <p(s, t) ^ (s + l)<p{X; s + 1, t - 1) 
or 

{ <p(s, t)^(t+ l)<p(s - 1, t + 1; Y) or <p(s, t) ^ (t + s - 1, t + 1) 



(7.8) { 



Proposition 7.7. Let p±, . . . ,pp be a family V of polynomials in 2k 2 commuting vari- 
ables. Then the SVR algorithm coupled with Alg.2 will yield an nc representation of 
p(X, Y) ofV with p in NC( 7 .$) if and only if the given family V admits an nc repre- 
sentation in this set. 

Proof. The proof is similar to the proof of Proposition 7.6, except that Alg.2 is used in 
place of Alg.l. □ 

7.2.2. Two letter monomials ^-equivalent to x s u x\ with s = t + 1. The result for this 
case is simpler than the case when s > t + 2 > 4 because the condition 

(7.9) s<p(s,t) ^ (t + l)<p(s - l,t + 1) 

required for Alg.l to work after an application of the SVR algorithm reduces to 

<p{t + l,t) ^ <p(t,t + l) 

when s = t + 1. However, we always require that <p(t+l,t) ^ <p(t,t + l) so that we can 
successfully partition the diagonal variables. Therefore our conditions to insure that 
we can successfully partition the dyslexic diagonal pairs imply that Alg.l will always 
be successful. 

Let NC (j. io) denote the class of all nc polynomials p(X, Y) of degree d such that 
for some integer t with 



(7.10) 



t>2, <p(t + l,t)^0, (p(t+ l,t) ^ <p(t,t+ 1) and 
<p(t + 1, t) ^ the coefficient of (XYfX or the coefficient of (YXfY. 
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Proposition 7.8. Let px, . . . ,pk 2 be a family V of polynomials in 2k 2 commuting vari- 
ables. Then the SVR algorithm coupled with Alg.l will yield an nc representation of 
p(X, Y) ofV with p in NC( 7 .xo) if and only if the given family V admits an nc repre- 
sentation in this set. 

Proof. The assumption that cp(t + 1, t) ^ 0, cp(t + 1, t) ^ fx, (p(t + 1, t) ^ f2, and 
(pit + 7^ (p(t, t + 1) allow Lemma 2.17 to determine the partitioned diagonal 
pairs. The SVR Algorithm can then be employed. The above discussion implies that if 
we differentiate with respect to Y in the SVR Algorithm that the DiagParl Algorithm 
will successfully partition the remaining variables between X and Y. Furthermore, 
the ParPosX Algorithm will position X and once X is determined, Algorithm PosPol 
positions the polynomials in the reduced family, which positions the polynomials in 
p(X, Y). Once the polynomials are positioned, then Algorithm PosY positions Y and 
Algorithm NcCoef will successfully determine p(X, Y). □ 

Remark 7.9. We do not write out an analogous result based on Alg.2 because in order 
to partition the diagonal elements in Alg.2 we require that tp(t + 1, t) ^ <p(t, t + 1), 
which insures that Alg.l will be successful. 

7.2.3. Two letter monomials \>-equivalent to x s u x l v with s = t. The preceding two 
cases dealt with families of polynomials containing two-letter monomials > -equivalent 
to x\x l j with s > t > 2. In those cases Lemma 2.16 or Lemma 2.17 was applied first 
to determine the dyslexic diagonal pairs. Then Lemma 2.19 was applied to partition 
the diagonal elements. However, if the only two-letter monomials are > -equivalent 
to xfx^ with s = t, then Lemma 2.19 is not applicable. This section is devoted to 
developing an algorithm for partitioning the diagonal entries and determining an nc 
representation in this particular case. The main result of this subsection is Theorem 
7.11. It is convenient, however, to first establish a preliminary lemma: 

Lemma 7.10. Let px, ■ ■ ■ ,Pk 2 be a family V of polynomials in 2k 2 commuting vari- 
ables xx , ■ • • , x 2 k 2 with an nc representation p(X, Y). Let Xi and Xj be a dyslexic diag- 
onal pair and let sx, tx and s 2 , t 2 be pairs of positive integers such that Sx+tx = s 2 + t 2 
and either Sx = s 2 or sx + 1 = s 2 and suppose that p(X,Y) satisfies the following 
conditions when s + t = sx + tx: 



(1) Ifs > 0, then (p(s, t + 1; X) ^ if and only if<p(s, t + 1; Y) ^ 0. 

(2) Ifs > 0, then (p(X; s, t + 1) ^ if and only ifip(Y; s, t + 1) ^ 0. 

(3) ip{0, s + t + 1) = <p(s + t + 1, 0) = 0. 



Suppose further that and fix^x^Xy appear in some polynomial p r E V 

where q^0 and (3 ^ and let t and m be the largest integers such that 
and x1 2 ~ m x t2+m x v appear in p r . Then 




and 
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Proof. There are two steps: 

1. £ > m =>• x u E L(xi): If i > m and x u E 1 L(xi), then x u must belong to L(xj), 
i.e., either x u E R(xj) or x u E C(xj). But if x u E R(xj) and Xi and Xj are in the ss 
position of X and Y, respectively, then x u is in the st position of Y and x v is in the st 
position of X (as axl l x l ^x u and ftxfx^Xy are in the same polynomial p r ). Therefore, 
ip(si — £, t\ + £ + 1; Y) ^ 0. If Si = £, then this contradicts (3) and therefore is not 
a viable possibility. If si > I, then ip(si — £, t\ + £ + 1; Y) ^ and, (1) implies that 
5 = cp( Sl -£,t 1 + £ + l;X)^0. Thus, 5xf~ {l+l) x t l +l+1 x v belongs to p r , i.e., 



si = s 2 =>- <5x* 2 ( - £+1 ^x' 2+ ^ +1 x„ belongs to p. 



j 



Sl = S2 - l ==> <fof ~ ( * +2) z* 2+ *+ 2 x„ belongs to p r 

However, the definition of m implies that £ + 1 < m in the first case and £ + 2 < m 
in the second, both of which clearly contradict the assumption that I > m. Therefore, 

x u £ R(xj). 

A similar argument based on (2) serves to prove that x u ^ C(xj). Therefore, x u E 
claimed. This completes the proof of 1. 

2. m > I ==>- x v E L(xi): If m > I and x v E~ L(xi), then x v E L(xj). Suppose 
that in fact x v E R(xj). Then (p(s 2 — m,t 2 + m + 1;Y) ^ 0. If s 2 = then 
this contradicts (3); if s 2 > m, then </?(s 2 — m,t 2 + m + l;Y) ^ and, by (1), 
7 = V 9 ( s 2 - m, t 2 + m + 1; X) ^ and hence that the monomial 7x- 2 " m_1 x* 2+m+1 x M 
is in this polynomial, i.e., 

si — (m+l) t-i+m-i-l , -, 

Si = s 2 =>- v x / x « belongs top r 
Si = s 2 - 1 7a;- 1 ~ m x* 1+m x M belongs to p r 

Therefore, the definition of £ implies that m + 1 < £ in the first case and m < £ 
in the second, which clearly contradicts the assumption that £ < m. Consequently, 
x v E L{x>i} as claimed. A similar argument rules out the case x v E C(xj). 

□ 

Let NC( 7 ii) denote the class of nc polynomials p(X, Y) of degree d > 4 such that 
(7.11) 

r > 2, cp(r,r)^0, <p(2r, 0) = cp(0, 2r) = and 
if s > 0, t > and s + t = 2r, then t : X) ^ <p(s, t:Y)^0, 
if s > 0, t > and s + t = 2r, then y?(X; s, t) ^ ^(K; s, t) ^ 0. 

Theorem 7.11. Let pi, ■ ■ ■ ,p k 2 be a family V of polynomials in 2k 2 commuting vari- 
ables xx, ■ ■ ■ , x 2 fc2 of degree d > 4. Then Lemma 7.10, the SVR Algorithm and either 
Alg.l or Alg.2 will yield an nc representation p(X, Y) with p in the class NCijxi) if 
and only if the given family admits an nc representation in this class. 
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Proof. Under the given assumptions Lemma 2.8 may be applied to obtain the dyslexic 
diagonal pairs 

\ x ix j x ji }>■••! { x ik J X jk } • 

To ease the notation, we assume that x is and Xj s are in the ss position for s = 1, . . . , k 
and that one of these pairs is partitioned, i.e., for some fixed choice of s, Xi a G X and 

x j3 g Y. 

The objective is to determine L(xi a ) and L(xj s ) for each dyslexic diagonal pair. 
Then for s fixed, there must exist at least k — 1 other dyslexic diagonal variables 
{x h , x is-l , x ia+1 , • • • , x ik } which have the property that 

(7.12) L(x it ) n L(x is ) ^ for t ^ s. 

Equation (7.12) implies that the variables {x ir , ■ ■ ■ , Xi a _ x , Xi s+1 , ■ ■ ■ ,x ik } lie on the 
diagonal of X with x is and that the union of the L(x it ) contains all of the variables in 
X. Therefore, this process serves to partition the variables between X and Y. Once 
this is done, the SVR Algorithm and the final steps in Alg.l or Alg.2 beginning with 
Algorithm ParPosX will determine an nc representation for V. 

Let p st denote the polynomial in the st position with s ^ tin the array corresponding 
to p(X, Y). Then p st , will be of the form 

(7.13) s , 

+ x\x r ja (gx u + ex v ) + x r it x r jt (hx u + fx v ) H . 

Given that cp(r, r) ^ 0, we must have that cp(r, r; X) ^ and cp(r, r; Y) ^ 0. and that 
one of the following cases must hold: 

1. a^0,e^0 and c = g = 0, 

2. a = 0, e = 0, and c ^ 0, g ^ 0, 

3. a ^ 0, e ^ and c + 0, g ^ 0. 

Similarly, the assumption that <p(r, r) ^ implies that y?(X; r, r) ^ and <fi(Y; r, r) ^ 
implies that one the following must also hold: 

4. b ^ 0, / ^ and d = h = 0, 

5. 6 = 0, / = 0, and d ^ 0, h ^ 0, 

6. 6 ^ 0, / ^ and d ^ 0, /i ^ 0. 

If case 1 holds, we have that a = ip(r, r; X) and e = ip(r, r; Y), and consequently, that 
x u G R(xi s ) and x„ G R(xj s ). If case 2 holds we may conclude that c = <^(r, r; X) 
and g = (f(r,r; Y) and that G R(x is ) and x u G R(xj s ). Finally, if case 3 holds, 
we must resort to more subtle measures to partition x u and x v . Here we can apply 
Lemma 7.10 with s± = s 2 = r — 1 and t x = t 2 = r to partition x u and x v between 
R(x is ) and R(xj 3 ). Therefore, by ranging over t in the polynomials p st for a fixed s, 
we will be able to entirely determine both R(x is ) and R(xjJ using the above analysis. 

We may similarly analyze cases 4-6 and range over s in the polynomials p st for a 
fixed t to determine C(xi t ) and R(xj t ). Thus, by ranging over all polynomials in V, 
we will be able to determine L(xi a ) and L(xj 3 ) for each dyslexic pair x is and Xj s with 
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1 < s < k. By determining which of these sets have a nontrivial intersection, we will 
obtain a partitioning of the variables x\, ■ ■ ■ , x 2 fc 2 between X and Y. 

□ 

Remark 7.12. The indicated terms in p st in the first part of the preceding proof can 
also be grouped as 

Pst = Xuiax^x^ + bx r r l x r jt + gx\x r r x + hx^r 1 ) 

+x v {cxZ 1 x r 3a + dx r r x xl + exiy- 1 + fx^; 1 ) + ■■■. 

7.2.4. Two letter monomials ^-equivalent to x s u x v . The next result presents another 
way of determining an nc polynomial representation for families of polynomials con- 
taining two letter monomials of the form ax s u x v with s > 2 and a e E \ {0}. 
Let C7(7.i4) denote the class of nc polynomials of degree d > 4 such that 

<p{d-l,l)?<p{l,d-l), 
cp(d - 1, 1) ^ <p(d, 0) if 1,1)^0, 
(7.14) l(p(l,d-l)^(p(0,d) if (p(l,d-l)^0, 

V (Y;d- 1,1)^0, (p(d-l,l;Y)^0, 
y?(X;l,d-l) ^0, y>(l,d- 1;X) ^ 0. 

Proposition 7.13. Let pi, . . . ,Pk 2 be a family V of polynomials in 2k 2 commuting 
variables. Then the SVR algorithm will yield an nc representation of p(X, Y)ofV 
with p in iVC(7.i4) if and only if the given family V admits an nc representation in this 
set. 

Proof. Suppose first that (p(d — 1, 1) ^ 0. The assumptions in the last two lines of 

(7.14) guarantee that k 2 — k polynomials will contain four or more terms > -equivalent 
to x d ~ x x v with x v 7^ x u , whereas the assumption that cp(d—l, 1) ^ cp(l, d—1) guaran- 
tees that exactly k polynomials p it , . . . , p ik will contain either one or two monomials 
t> -equivalent to x d l ~ 1 x v with x u ^ x v . This allows us to identify them as diagonal en- 
tries and to partition them between X and Y. Therefore, we apply the SVR Algorithm 
to V by differentiating d — 1 times to obtain the family {T^~ l p\, ■ ■ ■ , T^~ l p k 2} with 
nc representation 

D d x j k p(X, Y) = (d- l)\<p(d, 0)X + (d- l)\<p(d - 1, 1)Y. 

Since (p(d, 0) ^ <p(d — 1, 1) and cp(d — 1, 1) ^ by assumption, we can partition the 
remaining k 2 — 2k variables between X and Y, 

The construction of an nc polynomial representation can now be completed by in- 
voking the algorithms ParPosX,PosPol, PosY and Algorithm NcCoef. The remaining 
case follows similarly. □ 

7.3. Summary of Two-letter Algorithms. We now summarize our results based on 
the analysis of two-letter monomials. The methods of this section bear on p(X, Y) of 
the form 

(7.15) p(X, Y) = d^XYf + d 2 (YXY + d 3 (XYfX + d A {YX) l Y + q(X, Y), 



NON-COMMUTATIVE REPRESENTATIONS 61 

where q(X, Y) is an nc polynomial containing no multiples of the first four monomials 
in (7.15). The effectiveness of the procedures is summarized by: 

Theorem 7.14. Let pi, . . . ,pi~2 be a family V of polynomials in 2k 2 commuting vari- 
ables Xi, . . . , x 2 k 2 and let Q denote the set ofnc polynomials p(X, Y) of degree d > 1 
that satisfy the properties in at least one of the following three lists: 

(1) For some s, t G N, s > t + 2 > 4, ip(s,t) ^ 0, and ip(s,t) ^ ip(t,s). 
Additionally, assume that p(X, Y) satisfies one of the following conditions: 

(1) t<p(s,t) ^ (s + iMs + M-i), 

(2) s<p(s,t) + (t+%(s-l,t + l), 

(3) <p(s,t) ^<f(s + l,t-l;X) or <p(s,t) ^ <p(X; s + l,t - 1), 

(4) <p(s,t) ^ <p(s- l,t + 1;Y) or <p(s,t) ^ ip{Y; s - l,t + 1). 

(2) For some t > 2, ip(t + 1, t) ^ and additionally assume that p(X, Y) satisfies 
the following properties: 

(1) <p(t + l,t) ^d 3 

(2) tp(t + l,t)^d4, 

(3) <p(t + l,t)^<p(t,t+l) 

(3) For some r > 2, (p(r, r) ^ and additionally assume that the nc representation 
also satisfies the following properties: 

(1) <^(0,2r)=^(2r,0) = 

(2) (f{r,r) ^ di 

(3) (f{r,r) ^ d 2 

(4) cp(s,t;X) ^ cp(s,t;Y) ^ for all s,t such that s > 0,t > 0, 
s + t = 2r. 

(5) (p(X;s,t) ^ (p(Y;s,t) ^ for all s,t such that s > 0,t > 0, 
s + t = 2r. 

Then the two letter algorithms developed in Section 7.2 determine an nc representation 
p(X, Y) for V in Q if and only ifV has a representation in the class Q. 

Proof. Conditions (1), (2) and (3) are exactly what was needed to make the algorithms 
described §7.2.1, §7. 2. 2, and §7.2.3 effective. □ 

We can now supply a proof that our algorithms work for a collection of families of 
commutative polynomials that is in direct correspondence with a generic subset of the 
space of nc polynomials. 

7.3.1. Proof of Theorem 1.2. Suppose that U is a subspace of the space W of nc poly- 
nomials of degree d > 4. Then U must contain an nc monomial of one of the following 
forms as a basis element: 

m a% p{X, Y) with a = (ax, • - • , a n ), (3 = (ft, ■ ■ • , ft) and |a| + |ft = d, 

where ax, . . . , a n and ft, . . . , ft are positive integers, except that ax and ft are per- 
mitted to be equal to zero. Additionally, a and (3 must satisfy one of the following 
conditions: 

(1) |a| = d and |ft =0 or |a| = and |ft = d, 
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(2) \a\ > 0, |/3| > and either \a\ > \(3\ + 1, or \(3\ > \a\ + 1, 

(3) |a| > 0, \/3\ > and \a\ = 

(4) |a| > 0, \(3\ > and either \a\ = 1 or \(3\ = 1. 

If case (1) occurs, and p G U is a polynomial only in the single variable X or F, 
then the single letter algorithms will determine p and the set 

5i = (JV(7 ( 4.25)U ^67(4.27) nW 

will be open and dense in W since the indicated constraints are inequalities. 
Similarly, in case (2) the set 

s 2 = {NC (7 . 7) uNC {7 . w) r\u 

will be dense in U 

If case (3) occurs and U. contains either the term X d or Y d as a basis element, then 
case (1) applies. If hi does not contain this term as a basis element, then the set 

s 3 = NC {r . u) nu 

is open and dense in hi. 

If case (4) occurs, then the set 

s A = nc (7A I) n u 

is an open dense set in U given that the defining constraints are inequality constraints. 

□ 

7.4. Uniqueness results for two-letter algorithms. The family Q defined in The- 
orem 7.14 provides us with a large collection of nc polynomials for which our two 
letter algorithms will be successful. We now investigate the uniqueness properties of 
families with a representation in Q. 

Theorem 7.15. Suppose that V is a family of k 2 polynomials in 2k 2 commuting vari- 
ables that admits two nc representations p(X, Y) and p(X, Y) in the family Q. Then 
the matrices X, Y and X, Y are permutation equivalent (as defined in (1.13)). 

Proof. Suppose that (1) of (1.13)) is satisfied and let diag{X} denote the diagonal 
entries of the matrix X. The assumption that both p and p are in Q implies that either 

(a) diag{X} = diag{X} and diag{F} = diag{F} or 

(b) diag{X} = diag{F} and diag{F} = diag{X}. 

In case (a), if rup(s,n) ^ (s + l)y?(s + l,n — 1), then Remark 7.3 implies that 
Dp Ik p(X, Y) (resp., D~ J p{X, Y)) is an nc representation of the family TyPx, ■ ■ ■ , TyPk 2 
(resp., T~pi, ■ ■ ■ , T~p k 2). Furthermore, given that case (a) holds, the diagonals of Y 

and Y are the same, which implies that the families Tpp 1 , ■ ■ ■ , Typ k i and T~pi, ■ ■ ■ , T~p k 2 
are the same. 

Therefore Dtj p(X, Y) and D s ~ p(X, Y) are both nc representations of the family 
Typi, ■ ■ ■ , T Y p k 2. Moreover, given that p,p G Q, the nc polynomials Dp Ik p(X, Y) 
and D s ~ p(X, Y) are both in the set W that is defined in Theorem 1.4. Therefore, we 
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may apply Theorem 1 .5 to conclude that there exists a permutation matrix II such that 
either X = U T XU and J = Il T YIT or X = U T X T U and Y = U T Y T U. 

In case (b), X = U T YU and Y = U T XU or X = U T Y T U and Y = IT T X T IL The 
other three cases in (1.13) are handled in much the same way. □ 
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